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Fourth Order Theories of Gravity have recently attracted a lot of interest as candidates to explain 
the observed cosmic acceleration, the flatness of the rotation curves of spiral galaxies, the large scale 
structure and other relevant astrophysical phenomena. This means that the "Dark Side" issue of the 
Universe could be completely reversed considering dark matter and dark energy as "shortcomings" 
of General Relativity in its simplest formulation (a linear theory in the Ricci scalar R, minimally 
coupled to the standard perfect fluid matter) and claiming for the "correct" theory of Gravity as that 
derived by matching the largest number of observational data, without imposing any theory a priori. 
^ ■ As a working hypothesis, accelerating behavior of cosmic fluid, large scale structure, potential of 

| galaxy clusters, rotation curves of spiral galaxies could be reproduced by means of extending General 

Relativity to generic actions containing higher order and non-minimally coupled terms in curvature 
invariants. In other words, gravity could act in different ways at different scales and the above 
"shortcomings" could be due to the incorrect extrapolations of the Einstein theory, actually tested 
at short scales and low energy regimes. Very likely, what we call "dark matter and "dark energy" 
could be nothing else but signals of the breakdown of General Relativity at large scales. Then, it is a 
crucial point testing these Extended Theories in the weak field limit. In this sense, comparing these 
theories to General Relativity could be a fundamental step to retain or rule out them. In this review 
paper, after a survey of what is intended for Extended Theories of Gravity in the so called metric 
approach, we extensively discuss their Newtonian and the post-Newtonian limits pointing out, in 
details, their resemblances and differences with respect to General Relativity. Particular emphasis 
is placed on the exact solutions and methods used to obtain them. Finally,it is clearly shown that 
General Relativity results, in the Solar System context, are easily recovered since Einstein theory is 
a particular case of this extended approach. This is a crucial point against several wrong results in 
literature stating that these theories (e.g. /(_R)-gravity) are not viable at local scales. 
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; I. INTRODUCTION 

o . 

In recent years, the effort to give a physical explanation to the today observed cosmic acceleration [l| has attracted 
a good amount of interest in Fourth Order Gravity (FOG) considered as a viable mechanism to explain the cosmic 
acceleration by extending the geometric sector of field equations without the introduction of dark matter and dark 
energy. At fundamental level, several efforts have been aimed towards the unification of gravity with the other 
interactions of physics, like Electromagnetism, assuming GR as the only fundamental theory capable of explaining 
the gravitational interaction. The failure of such attempts led to the common belief that GR had to be revised in 
the ultraviolet limit in order to address issues like quantization and renormalization. These are only some aspects 
of the several physical and mathematical motivations to enlarge GR to more general approaches. For comprehensive 
reviews of the problem, see @-0|- 

Other issues come from astrophysics. For example, the observed Pioneer anomaly problem [S] can be framed into 
the same approach @ and then, apart the cosmology and quantum field theory, a systematic analysis of such theories 
urges at small, medium and large scales. In particular, a delicate point is to address the weak field limit of any 
extended theory of gravity since two main issues are extremely relevant: i) preserving the results of GR al local scales 
since they well fit Solar System experiments and observations; ii) enclosing in a self-consistent and comprehensive 
picture phenomena as anomalous acceleration or dark matter at Galactic scales. 

It is straightforward to extend GR to theories with additional geometric degrees of freedom and several recent 
proposals focused on the old idea of modifying the gravitational Lagrangian in a purely metric framework, leading 
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to fourth-order and higher-order field equations [10l-ll3]. Such an approach has become a sort of paradigm in the 
study of gravitational interaction consisting, essentially, in adding higher order curvature invariants and minimally 
or non-minimally coupled scalar fields into dynamics which come out from the effective action of some unification or 
quantum gravity theory. 

The idea to extend Einstein's theory of gravitation is fruitful and economic also with respect to several attempts 
which try to solve problems by adding new and, most of times, unjustified ingredients in order to give self-consistent 
pictures of dynamics. The today observed accelerated expansion of the Hubble flow and the missing matter at 
astrophysical scales are primarily enclosed in these considerations. Both the issues could be solved by changing 
the gravitational sector, i.e. the l.h.s. of field equations. The philosophy is alternative to add new cosmic fluids 
(new components in the r.h.s. of field equations) which should give rise to clustered structures (dark matter) or 
to accelerated dynamics (dark energy) thanks to exotic equations of state. In particular, relaxing the hypothesis 
that gravitational Lagrangian has to be a linear function of the Ricci curvature scalar R, like in the Hilbert-Einstein 
formulation, one can take into account an effective action where the gravitational Lagrangian includes other scalar 
invariants. 

Due to the increased complexity of the field equations, the main body of theoretical works dealt with the effort to 
achieve some formally equivalent theories which could be handled in a simpler way. In this sense, a reduction of the 
differential order of the field equations can be achieved by considering metric and connection as independent objects 
in the so called Palatini approach [Til . [l5j and energy conditions have to be carefully discussed [l6j . 

In addition, other authors exploited the formal relations to Scalar- Tensor theories to make some statements about 
the weak field regime [l8j], which was already worked out for scalar-tensor theories [l9j]. Also a post-Newtonian 
parameterization with metric approach in the Jordan Frame has been considered (20j . 

In this review paper, we want to address the general problem of the weak field limit for theories of gravity where 
higher order curvature invariants are present. In particular, we deal with theories where Riemann tensor, Ricci tensor, 
and Ricci scalar are considered in the effective action. We deduce the field equations, discuss the weak field limit and 
derive the weak field potentials with corrections to the Newtonian potential. The plan of the paper is the following: 

In seciion[Tll we provide a short summary of GR and discuss its extension to FOG. Besides we take into account the 
conformal transformations showing how such theories can be discussed under the standard of scalar-tensor gravity. 

In section ULTl we give generalities on spherically symmetric, Birkhoff theorem, Eddington parameters and deviations 
from GR solutions. A general perturbation approach for /(i?)-gravity is discussed starting from GR. 

Section [IV] is devoted to the technical development of field equations with respect to Newtonian and post-Newtonian 
approaches. 

In section [V] we take into account the debate on the analogy between f(R)- gravity and Scalar- Tensor gravity. 
We show that /(i?)-models are dynamically equivalent to O'Hanlon models which is a special case of Scalar- Tensor 
gravity characterized by a self-interaction potential without a kinetic term. We show that comparing the results of 
this theory with GR, in the weak field limit, could lead to misleading conclusions. 

In section IVI1 we analyze the Newtonian limit of J(i?)-Gravity. We are going to focus exclusively on the weak 
field limit within the metric approach. By using the development for a generic analytic function f(R), we find the 
solution in the vacuum with standard coordinates. Besides, we show that the Birkhoff theorem is not a general result 
for /(i?)-gravity since time-dependent evolution for spherically symmetric solutions can be achieved according to the 
order of perturbation. In other words, solutions could not be stable as requested by the Birkhoff theorem. 

In section IVII1 we report the Newtonian and post-Newtonian limit of field equations by using the Green function 
method for /(i?)-gravity and the formal solutions in the harmonic gauge. A general discussion about the mathematical 
properties of equations, the solutions and their deviations with respect to Gauss and Birkhoff theorem, and the 
Minkowskian behavior of metric tensor are reported. We derive the general solutions (at Newtonian and post- 
Newtonian levels) when an uniform massive spherical source is considered. The point-like source limit of the Newtonian 
solution is discussed and the compatibility of /(-R)-gravity with respect to GR is shown. 

In section IVIII1 we discuss the Newtonian limit but for a quadratic gravity Lagrangian where quadratic curvature 
invariants are presents. Also in this case, we adopt the Green function method. We find the internal and external 
potential generated by an extended spherically symmetric matter source. A particularly detailed discussion about the 
field equations and their solutions is provided for a set of values of the arbitrary constants present in the model. 

In section ILXl the Newtonian limit of the most general FOG is discussed with no gauge condition. The most general 
theory with fourth order differential equations is obtained by generalizing the /(-R)-models in the action with a generic 
function containing the other two curvature invariants: Ricci square (R a pR a P) and Riemann square {R a p 1 sR a ^ lS )- 
Considering the Gauss - Bonnet invariant, it is possible to show that only two of these invariants are really necessary. 

In section [Xj we draw the conclusions and discuss the possible applications of the results. 
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II. GENERAL RELATIVITY AND ITS EXTENSIONS 

Any relativistic theory of gravity has to match some minimal requirements to address gravitational dynamics. First 
of all, it has to explain issues coming from Celestial Mechanics as the planetary orbits, the potential of self-gravitating 
systems, the Solar System stability. 

This means that it has to reproduce the Newtonian dynamics in the weak held limit and then it has to pass the 
Solar System experiments which are all well founded and constitute the test bed of GR [2l[ . 

Besides, any theory of gravity has to be consistent with stellar structures and galactic dynamics considering the 
observed baryonic constituents (e.g. luminous components as stars, sub-luminous components as planets, dust and 
gas), radiation and Newtonian potential which is, by assumption, extrapolated to galactic scales. 

The third step is cosmology and large scale structure which means to reproduce, in a self-consistent way, the 
cosmological parameters as the expansion rate, the Hubble constant, the density parameter and the clustering of 
galaxies. Observations probe the standard baryonic matter, the radiation and an attractive overall interaction, acting 
at all scales and depending on distance. From a phenomenological point of view this is gravity. 

GR is the simplest theory which partially satisfies the above requirements [221 ]. It is based on the assumption that 
space and time are entangled into a single spacetime structure, which, in the limit of no gravitational forces, has to 
reproduce the Minkowski spacetime. Besides, the Universe is assumed to be a curved manifold and the curvature 
depends on mass-energy distribution [23j . In other words, the distribution of matter influences point by point the 
local curvature of the spacetime structure. 

Furthermore, GR is based on three first principles that are Relativity, Equivalence, and General Covariance(see 
0; S HI| for detailed discussions) . Another requirement is the Principle of Causality that means that each point of 
spacetime admits a notion of past, present and future. 

Let us also recall that the Newtonian theory, the weak held limit of GR, requires absolute concepts of space and 
time, that particles move in a preferred inertial frame following curved trajectories function of the sources (i.e., the 
"forces"). 

On these bases, GR postulates that gravitational forces have to be expressed by the curvature of a metric tensor 
field ds 2 = g a pdx a dx^ on a four-dimensional spacetime manifold, having the same signature of Minkowski metric, 

here assumed to be (H ). Curvature is locally determined by the distribution of the sources, that is, being the 

spacetime a continuum, it is possible to define a stress-energy tensor T^ v which is the source of the curvature. 

Once a metric g^ v is given, the inverse g^ v satisfies the condition 1 

g^gap = K (i) 

Its curvature is expressed by the Riemann tensor (curvature) 

where the comas are partial derivatives. The r^„ are the Christoffcl symbols given by 

r JL = t;.9 qct (5W + - 9hv,<t) (3) 
if the Levi-Civita connection is assumed. The contraction of the Riemann tensor ^ 

R^v — R ^ av — -T M „ iCT — ^ fie,i> + ^ fiu^ap ~ r^r^ (4) 
is the Ricci tensor and the scalar 

R = g^R aT = R\= g^T^ a - g^T^ ( + g^T^ - g^T^ (5) 

is called the scalar curvature of g^ v . The Riemann tensor ^ satisfies the so-called Bianchi identities and the 
contracted Bianchi identities, that is 



1 The Greek index runs between and 3; the Latin index between 1 and 3. 
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where the covariant derivative is A a ^--- S . 



V Q V" = □ = 



V ^ A a0...S = A ap...S^ +T a^ A <T/}...5 . 



. p/3 4a<r...5 . 



(6) 



r* /x A a ' 3 - ff and 



— is the d'Alembert operator with respect to the metric g^ v (see for the details [25[). 



Assuming that matter is given as a perfect fluid, that is 



T, 



(7) 



where u^u u defines a comoving observer (with the conditions g u Ut 2 — 1, = where x° — ct and we are assuming 
natural units with c = 1; p is the pressure and p the mass-energy density of the fluid, then the continuity equation 
requires T^ v to be covariantly constant, i.e. it has to satisfy the conservation law 



that are nothing else but contracted Bianchi identities. The GR field equations are then 



(8) 



G 



where 



X Tn, V 



R 



(9) 



(10) 



is the "Einstein tensor'''' of g M „. These equations are both variational and satisfy the conservation laws ((HJ since the 
following relation holds 



G p,a .„ = 



(11) 



as a byproduct of the above Bianchi identities [25|, l26[ . 

The Hilbert-Einstein Lagrangian that allows to obtain the field equations Q is the sum of an ordinary matter 
Lagrangian C m (minimally coupled) and of the Ricci scalar: 



£-he — \—g{R + XC r , 



(12) 



where \J—g denotes the square root of the value of the determinant of the metric g^ and the coupling constant is 
X = 8ttG. The action of GR is 



A GR = J d 4 xy/=g~{R + XC m ) 
From the action principle, we get the field equations (J9j) by the variation: 



SA GR = 8 J d i xy/=g(R + XC m ) = J d i x^g 
where T^„ is energy momentum tensor of matter 



R 

Rp,IS X T^ u 



(13) 



J d 4 xV^gg^8R^ = o (14) 
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The last term in (|14[) is a 4-divergence 

*Xy/=ggrSRp» = J d A x^—g[{-5g^),^ - D(g^Sg^)} (16) 

then we can neglect it and we get the field equation (j9|). For the variational calculus (Q~4|), the following relations can 
be used 

f = g^Sg^ + g^Sg^ = 

<V~S = -5 \^99apSg al3 

(17) 

<5i? = R a0 8g a P + g a P8R a p 
^ 5R af} = Sg p (a . l3}p - \U5g al3 - lg pa 5g p<7 . afi 

The choice by Hilbert and Einstein was completely arbitrary but it was the simplest one both from the mathematical 
and the physical viewpoint. As it was later clarified by Levi-Civita curvature is not a "purely metric notion" but, 
rather, a notion related to the "linear connection" to which "parallel transport" and "covariant derivation" refer [27| . 

This is the precursor idea of the " gauge theoretical framework" 28] , defined after the pioneering work by Cartan 
29]. 

It was later clarified that the principles of relativity, equivalence, covariance, and causality, just require that the 
spacetime structure has to be determined by either one or both of two fields, a Lorentzian metric g and a linear 
connection T, assumed to be torsionless for the sake of simplicity. 

The metric g fixes the causal structure of spacetime (the light cones) as well as its metric relations (clocks and 
rods); the connection T fixes the free-fall, i.e. the locally inertial observers. They have, of course, to satisfy a number 
of compatibility relations which amount to require that photons follow the null geodesies of T, so that T and g can 
be independent, a priori, but constrained, a posteriori, by some physical restrictions. These, however, do not impose 
that r has necessarily to be the Levi-Civita connection of g [3(| ■ 

This achievement justifies the fact, at least in principle, "Extended Theories of Gravity" can be considered starting 
from the same points by Einstein and Hilbert: they are theories in which gravitation is described by either a metric 
(the so-called "purely metric theories"), or by a linear connection (the so-called "purely affine theories") or by both 
fields (the so-called " metric-affine theories", also known as "first order formalism theories"). In these theories, the 
Lagrangian is a scalar density of the curvature invariants constructed out of both g and T. 

The choice (fl~2|) is by no means unique and it turns out that the Hilbert-Einstein Lagrangian is in fact the only choice 
that produces an invariant that is linear in second derivatives of the metric (or first derivatives of the connection). 
A Lagrangian that, unfortunately, is rather singular from the Hamiltonian viewpoint, in much than same way as 
Lagrangians, linear in canonical momenta, are rather singular in Classical Mechanics (see e.g. [3l|V 

A number of attempts to generalize GR (and unify it to Electromagnetism) along these lines were followed by 
Einstein himself and many others (Eddington, Weyl, Schrodinger, just to quote the main contributors; see, e.g., [32] ) 
but they were given up because of a number of difficulties related to the complicated structure of a non-linear theory 
(where by "non- linear" we mean here a theory that is based on non- linear invariants of the curvature tensor), and 
also because of the new understanding of physics that is currently based on four fundamental forces and requires a 
general "gauge framework" to be adopted (see [Hj]). 

Further curvature invariants or non-linear functions of them should be also considered, especially in view of the fact 
that they have to be included in both the semi-classical expansion of any quantum Lagrangian or in the low-energy 
limit of a string Lagrangian. 

Moreover, it is clear from the recent astrophysical observations that Einstein equations are no longer a good test 
for gravitation at all scales, that is at Solar System, galactic, extra-galactic and cosmic scales, unless one does not 
admit that the r.h.s. of Eqs.© contains some kind of exotic matter-energy density which are generically addressed 
as the " dark matter" and " dark energy" . 

From our point of view, the philosophy of Extended Theories of Gravity is much simpler (7| . Instead of changing 
the matter side of Einstein Equations ^ in order to fit the " missing matter-energy" content of the currently observed 
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Universe by adding new forms of matter and energy, we assume that it is more convenient to change the gravitational 
side of the equations, admitting corrections coming from non-linear Lagrangians. However such an approach should 
be explored. Of course, the Lagrangians should be conveniently tuned up on the basis of their best fit with all possible 
observational tests, at infrared scales (Solar, galactic, extragalactic and cosmic) and on the basis of consistency with 
fundamental theories at ultraviolet scales. 

Let us take into account an important class of Extended Theories of Gravity, the Fourth Order Gravity. It is the 
first straightforward generalization of GR and can be directly related to the ultraviolet (quantum gravity) and infrared 
(cosmology) issues of any final comprehensive theory of gravity. 

A. Fourth Order Gravity 

Let us start with the general class of Fourth Order Gravity (FOG) given by the action 



A* 



d X\J—g 



f(X,Y,Z) + XC r , 



(18) 



where / is an unspecified analytical function of curvature invariants X = R (Ricci scalar), Y = RapR" 13 (Ricci 
tensor square) and Z = R a p 1 gR a ^" fS (Riemann square). By varying the action (1181) in the metric approach and by 
using the properties (JTTJ) with following additional properties 



we have 



' SY = 8{R afs R afs ) = 2R IM a R au 5g^ + 2R^8R I1U 
5Z = 8{R aM5 R a ^ 5 ) = AR mPl R v ^8g^ + 2R a ^ 5 8R aM& 

8 R a 0^S = S(ga<jR a /3js) — R G 0~tS& 9aa + 9u<j8 R a fj lS 



8R a 



MS 



(/3;5)7 



Ph s] 



(19) 



8A f = 8 J d A xJ—g [f(X, Y, Z) + XC m ] 



d x-s/—g 



/ 



g^-XT^ )8g^ 



v fx + 2f Y R^)SR llv + 2f z R a ^ s 8R a 



M& 



J d 4 x^ j [f x R^ + 2f Y R«R av + ^f z R mM R v a ^ - f - 9lxv - X 
+fx[-(Sgn-^ - □(.9^%,)] + f Y R^[28g p {ix . v)p -USg^ - g" a 8 g„] 

+2f Z R a ^ l5 [8g a {i3-S) 1 - ^q(^ ;7 )5 + fig0{ T ,aS]] 



(20) 



~ J d A x^g-^f x R^ + 2frR«R av + 2f z R m&1 R v aM - f -g^ — X T^ v — f x .^ + g^Uf x 

-2[f Y R a (tl U a + U[f Y R^] + [f Y R aP \^g^ - m z R^J;a^Sg^ 

= J d 4 x^(H^ -XT^W W = 

where fx = dx> f Y = irr> fz = anc ^ the symbol ~ means that we neglected pure divergences. Then the field 
equations of FOG are 
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Hpv — fxRuv - -^9nv ~ fx-,nv + g^iuOfx + 2f Y R, J a Rav ~ 2\f Y R a ^]- v ) a + ^[IyR^v] + [fyRaff] '"^g^iis 
+2fzR l iaP^Rv al31 — ^[JzR^ a ^ v ];a/3 = X T^ v 

The trace of field equations (|2"Tj) is the following 



(21) 



H = g a ?H aP = f x R + 2f Y R a pR aP + 2f z R a ^ 5 R a ^ s - 2/ + D[3f x + f Y R] + 2[(f Y + 2f z )R a % a0 = XT (22) 

where T = T a a is the trace of energy-momentum tensor. Moreover the (|2"Tj) satisfies the condition H a>1 . a — 
XT atl . a = 0. In fact it is easy to check (in the case of /(-R)-Gravity) that 



H a ^, a = f[ a R a » + f'R a »., a - ^f'" - P a \ + P a a » = f"R a "R; a - f' a \ + f' a a M 
f"R a,J, R. a — f' a R a = f"R ail R; a - f"R' a R a " = 

where we used the properties G af * . a = and [V, V Q ]/' ;Q = - f a R a ^. 



(23) 



B. Conformal transformations 

The above theories can be easily viewed as scalar-tensor theories of gravity (for a comprehensive discussion see 0]). 
Let us now introduce conformal transformations to show that any scalar-tensor theory, in absence of ordinary matter, 
e.g. a perfect fluid, is conformally equivalent to an Einstein theory plus minimally coupled scalar fields. If standard 
matter is present, conformal transformations allow to transfer non-minimal coupling to the matter component [33.[35|. 
A general non-minimally coupled scalar-tensor theory is 



A 



ST 



d A Xyf^[F{(j))R + oj((f>)<p, a (P' a + V(<P) + XC r 



(24) 



where V((f>) and F(4>) are generic functions describing the self-interaction potential and the coupling of a scalar field 
4>, respectively. The Brans-Dicke theory of gravity is a particular case of the action (|2"4")l in which V(<p) = and 
lu(4>) — — iiL § s -. In fact we have 



A 



BD 



4>R - U>BD- 



XC r 



The variation of (|24p with respect to g^ and cj) gives the second-order field equations 



(25) 



F{4>)^ v + g^UFW = XT„ 



(26) 



2u((f>)a ( f>-u, ( b((f > )^ a ^ a - [F(4>)R + V(<I>)U = o 

3nF(4>)-F(<p)R- 2V{4>) - Lo{<p)(j), a <j)' a — XT 

I 2wfaP<£ + 3DF(4>) - \uj{4>) + u(4>)]4>. a ^ a - [F{<t>)R + V{4>))^ - F{<t>)R - 2V{4>) = XT 

The third equation in (|26[) is the trace of the field equation for g^ v and the last is a combination of the trace and of 
the field equation for <f). 

The conformal transformation on the metric g^ is 



g^v = A{x x )g lxv 



(27) 
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with A(x x ) > 0. A is the conformal factor. Obviously the transformation rule for the contravariant metric tensor is 
g^ v = A~ 1 g^ v . The mathematical quantities in the Einstein frame (EF) (quantities referred to g^) are linked to the 
ones in the Jordan Frame (JF) (quantities referred to g^ and action (|24|) ) as follows 



where 



' rj, = 1%, + + - ^9^ 

R^u = Rfiv - 20 ;M „ + 2(f> ;fl <p ;u - g^vOcf) - 2g^4> ;a (j)' a ' 
R = e-^^R - 6D0 - H;a<P a ) 

<t>;nv = 4>.,nv -^ a afi4>.a = (t>;nu - 1<t>;n<t>;v + 9tiv4>' a <t>;a 

In A 1 / 2 . The inverse relations are 



(28) 



(29) 



R a ^ v = R a ^ v + <^(M + <f>-,M - Wf^ + <t>-,M + 9^{4>~ a - <f> a <t>;p) - 9^~ a « - Tq <M 

R^v = Rfj,u + 2<j)- tlxv + 2(/)- <fl <p- u + g^ v 0<f> - 2g^4>' a (f). cr 
R = e 2 ^(R + 6D(j) - 60; ff ;(7 ) 

4>:p.v = 4>,p,v ~ r^^.CT = <f>- tl iv + 2<f>;n<f>;v ~ 9^^ ' 4>\a 

^ G M „ = G^v + 2(j).^ v + 2(j) ;fl (j). }U - 2<7^„ D(j> + g^fi"^ 

where □ and □ are the d'Alembert operators with respect to the metric <? M „ and g^ u . The transformation between 
the d'Alembert operators is □ = e 2tt >D - 2(f>' u d v . 

Under these transformations, the action (l24l) can be reformulated as follows 



A%jr = / <?Xy/-g 



AR+ VL{ip)ip, a ^ a + W{ip) + XL r , 



(30) 



in which R is the Ricci scalar relative to the metric g and A is a generic constant. The relations between the quantities 
in two frames are 



t»(<t>) _ 3 ( d\nF (0) 
F(4>) 2 



n(ip)df 2 = a 



d<p 2 



(31) 



r — 



A 



c, 



k F{4>)A{x x ) 1 = A 



The field equations for the new fields <7 M „ and ip are 
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< 2Q.(<p)n<p - Sl,< P ((p)<p; a (p ia - W tV (<p) = XC m>9 (32) 
R=-{ \Xf* + 2W(ip) + n{Lp)~g aT Lp, aV , T 

Therefore, any non-minimally coupled scalar-tensor theory, in absence of ordinary matter, is conformally equivalent 
to an Einstein theory, being the conformal transformation and the potential suitably defined by (|31l) . The converse is 
also true: for a given F ((/>), such that the relations (|31[) hold. we can transform a standard Einstein theory plus scalar 
field into a non-minimally coupled scalar-tensor theory. This means that, in principle, if we are able to solve the field 
equations in the framework of the Einstein theory in presence of a scalar field with a given potential, we should be 
able to get solutions for the scalar-tensor theories, assigned by the coupling F '((/>), via the conformal transformation 
(|2T|) with the constraints given by (|3"Tj) . Following the standard terminology, the "Einstein frame" is the framework of 
the Einstein theory with the minimal coupling and the "Jordan frame" is the framework of the non-minimally coupled 
theory [34], HH|. These considerations will be extremely useful for the discussion below. 

III. SPHERICAL SYMMETRY 

Spherical symmetry is the first step necessary to develop the Newtonian and the post-Newtonian limits of any 
alternative theory of gravity. Here, we will discuss the basic features related to such a symmetry starting from GR 
and then we develop a perturbation approach for spherically symmetric solutions of /(i?)-gravity. This machinery 
will be extremely useful to develop the FOG weak-field limit. 

A. Generalities on spherical symmetry 

Since we are interested to understand the modifications of GR predictions when one takes into account concentra- 
tions of matter in the space, it is fundamental to discuss properties of the metric g^. As a first step, we will study 
the gravitational potential generated by spherically symmetric distributions of matter. The most general spherically 
symmetric metric 2 can be written as 

ds 2 = gi (t, |x|) dt 2 + g 2 (t, |x|) dt x • dx + g 3 (t, |x|)(x • dx) 2 + g 4 (t, |x|)d|x| 2 (33) 

where gi are functions of the spatial distance |x| and of the time t. The set of coordinates is a; M = (t, x 1 , x 2 , a; 3 ). 
The scalar product is defined as x • dx — x 1 dx 1 + x 2 dx 2 + x 3 dx 3 . By the spherically symmetric form of (|33"|) , it is 
convenient to replace x with spherical coordinates r, 9, <fi defined as 

x 1 = r sin (9 cos 0, x 2 — r sin (9 sin </> , x 3 = rcos0 (34) 
The line element (f33|) then becomes 

ds 2 = gi(t,r)^ 2 +rg 2 (t,r)dtdr + r 2 g 3 (t,r)dr 2 + g 4 (t,r)(dr 2 +r 2 dQ) (35) 
where dfl = d9 2 + sin 2 Odtf 2 is the solid angle. We are free to reset our clocks by defining the time coordinate 

t = t' + e(t',r) (36) 



2 The metric is spherically symmetric if it depends only on x and dx only through the rotational invariants dx 2 , x • dx and x 2 . 
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with e(t',r) an arbitrary function of t' and r. This allows us to eliminate the off-diagonal element gt r in the metric 
(|55|) by setting 



de(t',r) rg 2 (t',r) 



dr 



(37) 



The metric ([55]) becomes 



ds 2 = gi(t',r) 



\ de(t',rY 


2 dt 12 + 


dt' 





dr 2 + g 4 (t',r)r 2 dn 



Agi{t',r) 

where introducing a new metric coefhcients gtt(t',r), g rr (t',r) and gnn(t',r), we can recast Eq. (I38|) as follows 



(38) 



ds 2 = g tt {t' ,r) dt' 2 - g rr (t' ,r)dr 2 - gnn(t' ,r)dfl 
Introducing a new radial coordinate (r') by considering the further transformation 



(39) 



r' = (const)e V 9 "« (t ,r > 
it is possible to recast Eq. ([39| into the isotropic form (isotropic coordinates) 



(40) 



ds 2 = g tt {t, r') dt 2 - gij (t', r^dx'dx 1 



(41) 



and then it is possible to choose gnn(t',r) 



(this condition allows to obtain the standard definition of the 



circumference with radius r") and to have the metric (|39[) in the standard form (standard coordinates) 



ds 2 = gtt{t' \ r ") dt 2 — g rr (t' : r")dr" 2 — r" 2 d£l 



(42) 



Obviously the functions g t t(t',r") and g rr (t',r") are not the same of (f3"9"| . If we suppose gij(t',r') = Y(t' ,r')5ij , it is 
worth noticing that it is possible to pass from (|4Tj) to (1421) by the coordinate transformation 



r' = r'(r") = {const)e-l ' dr " 



VY( r ") 



(43) 



We can, then, state that the expressions (|39|). (|41j) and (|42|) are equivalent to the metric (|33|) and we can consider 
them without loss of generality as the most general definitions of a spherically symmetric metric compatible with 
a pseudo - Riemannian manifold without torsion. The choice of the metric form is only a practical issue useful to 
develop calculations. 



B. The Birkhoff theorem in General Relativity 

The Birkhoff theorem is an important result of GR which essentially states that stationary solutions are also static 
and stable. Specifically, the theorem holds for spherically symmetric solutions. 
Let us start our considerations by rewriting the metric (|4"2"j) as follows 

ds 2 = e^' r) dt 2 - e^ r) dr 2 - r 2 dQ (44) 
where we redifined the radial coordinate. The only non- vanishing components of metric tensor g M „ are 



gtt = e y ■ ' , g rr 



u(t.r) 

-e w ■ ; , gee 



9<P<P 



2 • 2 

-r sin i 



(45) 
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with functions /i(t, r) and v{t,r) that are to be determined by solving the field equations in GR ((9j). Since g^ u is 
diagonal, it is easy to write all the non- vanishing inverse components: 



tt = e -„(t,r) > 



rr — u(t.r) 6 

.9 = -e PV \ 



-2 • -2 

-r sin i 



Furthermore, the determinant of the metric tensor is 



g = -e^+^-'Vsm 2 ! 



so the invariant volume element is 



(46) 



(47) 



\/—g dr d6 d(f> = r 2 e ' 2 sin 8 drd9 1 
The only non- vanishing components of Christoffel symbols ([3]) are 



(48) 



r' 
1 u 



£>(t,r) 
2 



pt ^' (t,r) u(t,r)-v(t,r 

rr 2 



-/i(t,r) 



= - sin0 cos ( 
and the field equations (|9]) become 



v'(t,r) -pr _ v'{t,r) v(t,r)- u(t,r) 



1 tr 



i/(t,r) 
2 



p0 _ p0 _ 
1 r6 — 1 r0 ~ 

rl =coti 



u(t,r) 



T r ,, = -re-^rt sin 2 6> 



(49) 



1 _ e -/i(*,r) 



i _ //fer) 



= XT* 



A(*i^) c -uft,r1 = %T r 



-fi(t,r) 



v'{t,r) 1 



/t(t,r) + AferMMO 



= AT a = AT* 



If we suppose a stress-energy tensor jJJ induced by a point-like source with mass M 



(50) 



T = p(x) 

where p = M5(x), we obtain the Schwarzschild solution in standard coordinates 3 



(51) 



ds 2 = 



1 - ^ 



dr" 2 

dt 2 - ^-^ _ r "2 dfi 



(52) 



3 5(x) is the 3-dimensional Dirac 5-function. 
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where r g = 2GM is the Schwarzschild radius. 

Metric (|52p determines completely the gravitational field in the vacuum generated by a spherical matter density 
distribution. Furthermore the Schwarzschild solution is valid also when we consider a moving source with a spherical 
distribution. The spatial metric is determined by expression of spatial distance element 



dr" 2 

dl 2 = ——T£ + r" 2 dn (53) 

r" 

We have to note that, while the length of circumference with "radius" r" is the usual one 2irr" , the distance between 
two points on the same radius is given by the integral 



dr" 

" " (54) 



r 

this means that the space is curved. Besides we note that gu < 1, then, by the relation between the time coordinate 
t and the proper time r (dr = i/gudt), we get the condition 



dr < dt 



(55) 



At infinity, the coordinate time coincides with the physical time. We can state that when we are at a finite distance 
from the the mass, there is a slowdown of the time with respect to the time measured at infinity. 

In presence of matter the situation is the following. From the first equation in the (|50|) . when r — > 0, fi(t, r) has to 
vanish as r 2 ; otherwise T t t could have a singular point in the origin. By formally integrating the equation with the 
condition fi(t, r)\ r — o = 0, one gets 



(i(t,r) 



In 



i_ * r T \f 2 df 

r Jo 



(56) 



Beside the point-like case, it is easy to demonstrate that the pro prieties (JMJ), (I55t and fi(t, r) + v(t, r) < hold also in 
more general matter distributions with spherical symmetry [25j. If the gravitational field is generated by a spherical 
body with "radius" £, we have T t t = outside the body (r > £) and we can write 



In 



I - — [ ' T\f 2 df 



(57) 



obtaining the analogous expression of (1521) in the matter 



ds 2 =, 



dr 



dr 



112 



1 



r g (r") 



r" 2 dn 



(58) 



where we have introduced the Schwarzschild radius related to the quantity of matter included in the sphere with 
radius r": 



r g {r") = X T\f 2 df 
Jo 

Obviously when the distance is bigger than the radius of the body, the metric ([55)) is equal to ([5 
If we consider the transformation (|43|) . which in the case of Schwarzschild solution is 



(59) 



2r" - r g + 2 v /r" 2 - r g r" 



(60) 



it is possible to obtain the Schwarzschild solution (|52p in isotropic coordinates 
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ds 2 = 


'1 - 


r 9 - 

ir' 


1 

dt 2 - 


In 


T 9 ' 




.1 + 


Is. 

At' - 




Ar 1 



(dr' 2 + r' 2 dn) 



(61) 



In both cases, solutions (|5"2"j) and (1551) obviously satisfy the trace of field equations: R = —X T. Since in the 
vacuum the trace of matter tensor is vanishing (except in the origin, where the trace is proportional to <5(x)) we can 
state that the Schwarzschild solution is "Ricci flat": R — 0. 

If we add in the Hilbert - Einstein Lagrangian ([12")) a term like (—2y/—g A), with A a generic constant, the field 
equations (|9]) are modified as follows 



G 



Agw = XT, 



(62) 



In the case of a point-like source, we find the Schwarzschild - de Sitter solution 



ds 2 = 



r" 3 



dt 2 - 



dr h 



+ f r" 2 



- r" 2 dn 



(63) 



The trace of (El is 



R : = 4A - X T (64) 

from which we note that this solution does not admit solution in vacuum, since also in absence of ordinary matter 
(T^v = 0) we have a non-vanishing scalar curvature. The contribution is given by the cosmological constant A. Also 
in this case, analogous considerations as in ([58]) hold. 

Finally let us consider as source a radial and static electric field E = Qx/|x| 3 . We know that the Lagrangian of 
electromagnetic field is —j-F a pF a P where F a p is the electromagnetic tensor. Then, the Hilbert - Einstein Lagrangian 
is 



Che = V^9(R - t F a pF a ?) (65) 

and the Einstein equation © becomes 

<V = ~ (g^F af3 F a P - AF m F a u ) (66) 
The solution for a spherically symmetric system is the Reissner - Nordstrom solution 



ds 2 = 



91' 

r ll2 



dt 2 - 



dr h 



1 



- r" 2 dtt 



(67) 



In all the above cases, the Birkhoff theorem holds: The metric tensor generated in vacuum by a matter density 
distribution with a spherical symmetry is time-independent. Also a time-dependent source with a spherical symmetry 
produces a static metric. The curvature of spacetime in the matter, at a distance r from the origin, is proportional 
only to the matter inside the sphere of radius r. This conclusion is compatible with the Gauss theorem of Classical 
Mechanics. 



C. The Schwarzschild solution and the Eddington parameterization 

The Schwarzschild radius r g in the solution (|6ip is a scale-length induced by theory and represents a natural 
parameter to study the gravitational interaction with respect to the generating source. At radial distance r' where 
r g /r' <C 1, solution (Hxf) becomes 
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ds 2 



dt 2 



dr' 



r ' 2 dn 



(68) 



that is the standard situation presented by any self-gravitating system far from its critical radius. 

Since we are interested to investigate the deviations, induced by FOG, from the behavior (l68l) . it is useful to intro- 
duce the method taking into account such deviations with respect to GR. The standard approach is the Parameterized- 
Post- Newtonian (PPN) expansion of the Schwarzschild metric (joTj) . Eddington parameterized deviations with respect 
to GR considering a Taylor expansion in term of r g /r' and assumed that, in Solar System, the limit r g /r' -C 1 holds 
[2lj . The resulting metric is 



ds 2 



1 n r 3,P( r 9 



dt 2 - 



1 + 7^ 
r 



dr 12 + r' 2 dn 



(69) 



where a, /? and 7 are dimensionless parameters (Eddington parameters) which parameterize deviations with respect 
to GR. The reason to carry out this expansion up to the order (r g /r') 2 in g tt and only to the order (r g /r') in gij 
is that, in applications to Celestial Mechanics, always appears multiplied by an extra factor v 2 ^ M/r'). It is 
evident that the standard GR solution for a spherically symmetric gravitational system in vacuum, is obtained for 
a = /3 = 7 = 1 giving again the "approximated" Schwarzschild solution (|68[) . Actually, the parameter a can be 
settled to the unity due to the mass definition of the system itself [21] . As a consequence, the expanded metric (f69|) 
can be recast in the form : 



ds 2 



1 _ H + /3 ~ 7 

r" 2 



dt 2 



1 + 7^ 



dr' 



r" 2 dn 



(70) 



where we have restored the standard spherical coordinates by means of the transformation 



1 + ^. 

L ' ifi 



The parameters /?, 7 have a physical interpretation. The parameter 7 measures the amount of curvature of space 
generated by a body of mass M at radius r'. In fact, the spatial components of the Riemann curvature tensor are, at 
post-Newtonian order, 



independently of the gauge choice, where Nijki represents the geometric tensor properties (e.g. symmetries of the 
Riemann tensor and so on). On the other hand, the parameter j3 measures the amount of non-linearity (~ {rg/r') 2 ) 
in the gtt component of the metric. However, this statement is valid only in the standard post-Newtonian gauge. 

These considerations can be developed for any relativistic theory of gravity but, as we shall see below, the above 
results, valid in GR, cannot be simply extrapolated to any modified theory of gravity since misleading conclusions 
could be achieved. For example, in literature, there are some papers stating that f(R) 7^ R arc not viable models 
in Solar System since, by recasting the /(i?)-gravity as the O'Hanlon model, 7I/2 [47| in evident contrast with GR 
measurement giving 7 ~ 1. On the other hand, assuming f(R) = R 1+€ with e -> O-the standard GR would be 
hard to be recovered in the weak field limit and the Cauchy boundary conditions [HI, |49[ would be highly violated 
if the theory results switched from 7 = 1, for e = 0, to 7 = 1/2, for e =^ . As shown in (55|, the shortcoming can 
be solved by considering separately the weak field limits of f(R) and O'Hanlon theory which have to be confronted 
"after" the PPN-approximation. The conceptual reason is clear: the invariance properties of the theories are lost in 
the approximation process so the confront has to be performed assuming the same gauge and formalism. In other 
words, different theories have to be carefully confronted at general level and at approximated level without incautious 
extrapolations. This argument will be considered in detail below. 



D. Perturbing the spherically symmetric solutions of /(i?)-gravity 



Let us focus now on the differences of perturbing a spherically symmetric solution coming from a generic f(R)- 
gravity model with respect to the standard spherically symmetric solutions coming from the GR case f(R) — R. We 
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neglect, for the moment, the contributions of Ricci and Riemann square [36j which will be considered below. The 
search for solutions in FOG can be faced by means of a perturbation theory. 

The general approach is to start from analytical /(i?- function assuming that the background model slightly deviates 
from GR. Such a method can provide interesting results on astrophysical scales where spherically symmetric solutions, 
characterized by small values of the scalar curvature, can be taken into account. Below, we will consider the perturbing 
approach assuming that the background metric matches, at zero order, the GR solutions. 

The held equations (J^TJ) and (|22[) in /(_R)-gravity, if we develop also the covariant derivatives, become 



H = f'R-2f + n = XT 



(72) 



where the two quantities and H read 



9 



— f"'\ R,iiR,u — g^g aT R,aR.r 



n = 3/" 



9 aT M + 9 aT In \F I 9~ R,r + 9 aT R 



(73) 



Zf'"g aT R, a R,T 



where f l is the z-th derivative of / with respect to Ricci scalar R. 

In general, searching for solutions by a perturbation technique means to perturb the metric 



9n 



- (o) , (i) 



(74) 



This implies that the held equations (I72[) split, up to hrst order, in two parts. Besides, a perturbation on the metric 
acts also on the Ricci scalar R defined in Eq.© 



R ~ R^+R^ 

and then we can Taylor expand the analytic function f(R) about the background value of R, i.e.: 



(75) 



f = T r(fl (0) ) 

•I Z— / n n ! 



fi, = T r +2 (fl (0) ) 



R - ijO) 
R - ijW 
R - i?< 0) 



fJV(O) 



f (o) + f m R (i) + r^w 2 + H^L R m 



2 fiv(o) 

_1_ I 

6 



df 



f m + f »>(°) R (i) + fXll R w 2 = jm^ 



f" = E„ 



R - R<f>) 



f'"(0) , fiv(°)p(l) _ df" 



The zero order of held equations (|?2")) reads 



fl(0) d(0) _ 1(0) f (0) , o/(0) = T (o) 



where 



(76) 



(77) 



_ r (o) r^(0)^(0) _ g(o) g (o) P « R (o )R( oj I (78) 
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At first order, one has : 



with 



/ /(o) te4< )jR(1) 



+ Z?(0) _ lf(0) (1) 4.-1/(1) _ vji(l) 



1/(1) = _f''( )//?(l) WO)? D(l)_ r (l) P E>(0)_ o (0) 



+g^°RW + gW»°R<$ + g^ (ln^ ( >« + In- 



,p 



pa 



+9®" in ^ Jl<g> ) j - f m {R$nfp + B$R$ - g { $9 {0)p ° (R^RV 



~ g^gM^nM - g^g^R^^J) ~ f m R {1) \ R% 

-f v{0) R^!R^R^-g^g(°^R(^R^ 



(79) 



(80) 



(81) 



A part the analyticity, no hypothesis has been invoked here on the form of the function f(R). As a matter of fact, 
f(R) can be completely general. At this level, to solve the problem, it is required the zero order solution of (|77l) which, 
in general, could not be a GR solution. This problem can be overcome assuming the same order of perturbation on 
the f(R), that is: 



f(R) = R + T(R) (82) 
where JF(R) is a generic function of the Ricci scalar (this means to consider F{R) <C R) . Then we have 



/ = flW + R<» + , /' = 1 + J-' (0) , /" = J-" (0) , /"' = T 



m(0) 



(83) 



fill p(0) \ p(l) 

However the condition F(R) <C R has to imply the validity of the linear approximation Wr?oJ] — ^ !• This is 
demonstrated by assuming /' = 1 + T' and /" = T" . Immediately, we obtain that the condition is fulfilled for 



< 1 



For example, given a Lagrangian of the form f(R) = R + Eq. ([84l) becomes 



(84) 



2i? i? (1 



R(°)(R(V -Ro) 



< 1 



(85) 



while, for /(i?) = i? + ai? 2 , Eq.([84]) becomes 



1 + 2ai?(°) 



<C 1 



(86) 
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This means that the validity of the approximation strictly depends on the form of the models and the value of the 
parameters, in the previous cases, Rq and a. For the considerations below, we will assume that condition (|84[) always 
holds. 

Eqs.([77|) reduce to the GR form 



On the other hand, Eqs. (J79)) reduce to 



(87) 



_ _ i„Wp(0) _ i„(0) F (0) , 77/(0) ff (0) .07(1) _ -y/Wl) 



where 



(89) 



The new system of field equations is evidently simpler than the starting one and once the zero order solution is 
obtained, the solutions at the first order correction can be easily achieved. 
By assuming the time-independent metric (|42]l 



ds 2 = a(r)dt 2 - b(r)dr 2 - r 2 dVi 2 (90) 

it is straightforward to obtain new spherically symmetric solutions by substituting into the above perturbed field 
equations. In Tables HI and [TT1 a list of solutions, obtained with this perturbation method, is given considering different 
classes of /(i?)-models. Some remarks on these solutions are in order at this point. In the case of /(i?)-models which 
are evidently corrections to the Hilbert-Einstein Lagrangian as A + R + eR In R and R + eR n , with e <C 1, one obtains 
exact solutions for the gravitational potentials a(r) and b(r) related by a(r) = fe(r) -1 . The first order expansion is 
straightforward as in GR. If the functions a(r) and b(r) are not related, for f(R) = A + R + eR In R, the first order 
system is directly solved without any prescription on the perturbation functions x{r) and y(r) indicated in the Tables. 
This is not the case for the models f(R) = R + eR n since, in this case, one can obtain an explicit constraint on the 
perturbation function. In such a case, no corrections are found with respect to the standard solution. The models 
f(R) = R n and f(R) = ^^^) show similar behaviors. The case f(R) = R 2 is peculiar and it has to be dealt 
independently. For details, see j45| . 



IV. GENERAL REMARKS ON NEWTONIAN AND POST-NEWTONIAN APPROXIMATIONS 



At this point, it is worth discussing some general issues on the Newtonian and post-Newtonian limits. Basically 
there are some general features one has to take into account when approaching these limits, whatever the underlying 
theory of gravitation is. In fact here we are not interested in entering the theoretical discussion on how to formulate a 
mathematically well founded Newtonian and post-Newtonian limits of general relativistic field theories, nevertheless 
we point the interested reader to [37l - l43l | . In this section, we provide the explicit form of the various quantities needed 
to compute the approximations in the field equations in GR theory and any metric theory of gravity. We only mention 
that there is also a discussion on alternative ways to define the Newtonian and post-Newtonian limits of FOG in the 
recent literature, see for example In this work, the Newtonian and post-Newtonian limits are identified by the 
maximally symmetric solutions, which are not necessarily the Minkowski spacetime in /(i?)-gravity which could be 
singular. 

Let us start with some very general considerations on the Newtonian gravitating systems. If one takes into account 
a system of gravitationally interacting particles of mass M, the kinetic energy ^Mv 2 will be, roughly, of the same 
order of magnitude as the typical potential energy U = GM 2 /f, with M, f, and v the typical average values of 
masses, separations, and velocities of these particles. As a consequence: 
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/ - theory: 
spherical potentials: 
solutions: 
first order metric: 



A + R + eR In R 
a(r) = fe(r)- 1 = 1 + *l _ + Sx(r) 



solutions: 



x{r) 



k2_ , eA[ln(-2A)-l]r 2 
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o(r) = l-^-+Ji(r), 6(r) = — ^ + 6y(r) 



C (r) = (Ar 2 -6)<jfa + J 36rS ^_ 6) 



4<5(2A 2 r 4 - 15Ar 2 + 18)y(r) + r{36reA[log(-2A) - 1] 



+S(Ar 2 - 6)V (r)} 

fc 2 <5-6r 3 eA[ln(-2A)-l] 



/ - theory: 
spherical potentials: 
solutions: 
first order metric: 
solutions: 



R+eR n 
a(r) = bir)- 1 = 1 + ^ + Sx(r) 
x(r) = !f 

( P ) = 1 + J-fel i 6 ( r ) = l + 5 2M 

x(r) = fa + far, y(r) = fa 



/ - theory: 
first order metric: 

solutions: 



R/(Ro + R) 



a;(r) = 
y(r) = 



4e Ve 



lfl_ r 



fa - 2 ^ e a /f fa + far 



^ (6-Ro /2 +yifl r) 
5W 2 



fa 



2e >/« (\/6--Rn /2 r) , 
„3/2 «2 



TABLE I: A list of exact solutions obtained via the perturbation approach for some classes of /(i?)-models; fa are integration 
constants; the potentials a(r) and 6(r) are defined by the metric (|90[1 , 



w 2 — 91 

r 

(for instance, a test particle in a circular orbit of radius r about a central mass M has velocity v, given in Newtonian 
Mechanics, by the exact formula v 2 — GM/r.) 

The post- Newtonian approximation can be described as a method for obtaining the motion of the system to higher 
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/ - theory: 
spherical potentials: 

solutions: 



R n 



a(r) = fe(r) 



1 + 



fcl I Rpr 
r ~ 12 



+ 5x(r) 



first order metric: 



solutions: 



first order metric: 



solutions: 



"O'-Q In(r-rp) 
8 + 3«0>~o 



2 r 

= 2, i? / and a;(r) = + hzL_ + !sa f drr *)f dr _ 

with ro satisfying the condition 6fci + 8ro + Ror^ = 

,n > 2, System solved only whit Ro = and no prescriptions on x(r) 



"n = 2 y (r) = -^-^ + \rx'(r) + k 1 , 
with R(r) = SRo 

n 2 y(r) = -§ / drr 2 R(r) - 2®. + \rx\r) + ki 
with -R(r) whatever 

o(r) = l-2t+fa(r), 6(r) = ^ + Sy(r) 



r\ f \ rfel I r &2 



r g r 2 x(r)+2(rgT , 3— r 4 )x'(r) 
(3i' g -4r)(r g -r) 2 



n/2 whatever functions a;(r) , y(r) and 



TABLE II: A list of exact solutions obtained via the perturbation approach for some classes of /(ii)-models; ki are integration 
constants; the potentials a(r) and b(r) are defined by the metric (|90[l . 



approximations than the first order 4 with respect to the quantities GM / f and v 2 assumed small with respect to the 
squared light speed. This approximation is sometimes referred to as an expansion in inverse powers of the light speed. 

The typical values of the Newtonian gravitational potential $ are nowhere larger (in modulus) than 10~ 5 in the 
Solar System (in geometrized units, $ is dimensionless). On the other hand, planetary velocities satisfy the condition 
v 2 < — <£>, while the matter pressure p, experienced inside the Sun and the planets, is generally smaller than the matter 
gravitational energy density — p&, in other words 5 p/ p < — $. Furthermore one has to consider that even other forms 
of energy in the Solar System (compressional energy, radiation, thermal energy, etc.) have small intensities and the 
specific energy density IT (the ratio of the ener gy density to the rest-mass density) is related to U by IT < U (II 
is ~ 10~ 5 in the Sun and ~ 10~ 9 in the Earth [2l|). As matter of fact, one can consider that these quantities, as 
function of the velocity, give second order contributions : 



_$^ w 2 ^p/ j0 ^n-O(2) (92) 

Therefore, the velocity v gives 0(1) terms in the velocity expansions, U 2 is of order 0(4), Uv of 0(3), UU is of 0(4), 
and so on. Considering these approximations, one has 



4 This approximation coincides with the Newtonian Mechanics. 

5 Typical values of p/p are ~ 10" 5 in the Sun and ~ 10" 10 in the Earth GW. 
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d_ 

at 



v 



(93) 



and 



\d/dt\ 



0(1) 



(94) 



The paradigm of post-Newtonian limit is to start from a series-developed metric tensor (and all additional fields 
in the theory) with respect to the dimensionless quantity v (in natural units!). A system of moving bodies radiates 
gravitational waves and thus loses energy. This loss appears only at the fifth-order in the approximation of v. In the 
first four approximations, the energy of the system remains constant. From this, it follows that a system of gravitating 
bodies can be correctly described by a Lagrangian up to terms of order v , in the absence of electromagnetic fields. 
We thus find the equations of motion of the system in the next approximation after the Newtonian. 

To solve the problem, we have to start from the determination of the metric tensor g^ v in the same approximation 
of the gravitational field [45|, [46| . 

Considering that particles move along the geodesies 



d 2 x^ „ dx° dx T 
1_ t;m _ q 

ds 2 CTr ds ds 



these can be written in details as 



d 2 x l 
~di 2 ~ 



= -n 



dx" 



dt 



dx m dx n 
'~dt ~d~t 



r 



2r 



dx r ' 



dt 



2Ti 



dx m dx n 



dt dt 



dx l 
~dt 



(95) 



(96) 



In the Newtonian approximation, that is vanishingly small velocities and only first-order terms in the difference 
between g^ u and the Minkowski metric 77^, one obtains that the particle motion equations reduce to the standard 
result 



dt 2 



1 dg tt 

2 dx l 



The quantity 1 — g tt is of order GM/r, so that the Newtonian approximation gives 
is, to the order v 2 jr. As a consequence, if we would like to search for the post-Newtonian approximation, we need 



d 2 x i 
~dP 



(97) 



to the order GM/f 2 , that 



to compute 



dt 2 



to the order v A jr. Due to the Equivalence Principle and the differentiability of spacetime manifold, 



we expect that it should be possible to find out a coordinate system in which the metric tensor is nearly equal to the 
Minkowski one rj^ Vl the correction being expandable in powers of GM/f ~ v 2 . In other words, one has to consider 
the metric developed as follows : 



f g tt (t, x) ~ 1 + g^it, x) + g { t t\t, x) + 0(6) 



5tl (t,x)~ 5t ( l 3) (t,x)+0(5) 

I g y (i,x)~~^-+4 2) (i,x) + 0(4) 
where 5{j is the Kronecker delta, and for the controvariant form of g^ v , one has 



(98) 



' g u {t, x) ~ 1 + g (2)tt (t, x) + g^ tt (t, x) + G(6) 

3 "(t,x)~ 3 ( 3 )"(t,x)+0(5) 
[ ff «(t,x)~-^+ ff ( 2 )y(t,x)+0(4) 



(99) 
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The inverse of the metric tensor (IMl) is defined by ([TJ. The relations among the higher than first order terms turn 
out to be 



r ff ( 2 )«(t,x) = - 5 ( t 2 »(i,x) 



( 4 )"(i,x)= 5 ^(t,x)"- 5 W(t,x) 



(4), 



7 (8)« = ,(3) 



(100) 



In evaluating T^a we must take into account that the scale of distance and time, in our systems, are respectively set 
by f and f/v, thus the space and time derivatives should be regarded as being of order 



d 

dx l 



d v 
dt r 



Using the above approximations ([98]). ([99]) and (|100j) we have, from the definition ((3]), 



P(2) 1 _ 1 [ ( 2 ) _ ( 2 ) _ ( 2 ) 
1 jfc — 2 I Hjfc,i Hy,fc Hjfc,j 



rO)V - ifo (3) -o (3) -a (2) 

1 tj ~~ 2 \ ytj.z Hit, j Hjj,t 



The Ricci tensor components (|4]) arc 



r( 3 )'. = ±fa (3) W 3) - (2) 

2 I Hit j ~ Hjt,j yij,< 



ti 2 



r< 4 ) 1 - ±(V 4) W 2 V 2) -2o (3)N l r< 2 )' 

tt — 2|Htt,i i yim,ytt,m ^iiit,t I 1 t 



1 (2) 
29tt,i 



(101) 



(102) 



.(2) _ 1 (2) 
tt ~ 2^tt,rr 



d(4) _ 1 (4) ,1.(2) J2) 1 (2) (2) 1 (2) i(2) (2) _ 1 (2) (2) (3) 

^tt — 2ytt,rnm ~<~ 2 H mn i m i»tt,fi 2H m «Htt,m« "T" 2ilrnm,tt 4Htt,mHtt,m 4H«"«,™Htt,n Htm,tm 



1? 



(3) _ 1 (3) 
Hi 



2 

1 J 2 ) 



1„(3) 



1„(2) 



_ i W I i fl W 

mi 2Hmt,mz 1 2Hmm.tz 



7? 



(2) _ 1 (2) 



1 (2) . 1 J2) 



Sj 2^J,mm 2 ^im-mj 2Hj'Tn,7m 2@tt,ij 2Hmm,ij 

and the Ricci scalar (|5|) is 



'.9, 



(103) 



f p(2) - p( 2 ) _ r( 2 ) - „( 2 ) _ J 2 ) i „ 

■n, — Xl tt rt, m rn — ytt,mm <jnn,min T </ 
R(4) - I?(4) _ (3) d(2) _ (2) „(2) 



(2) 

mn.mn 



2ytt,mm ~ 2y mn > m ytt,n ~ 2H"»™Htt,mn ~ jSmm,8 4a 'tt ,mi 'tt ,m T 

_I ( 2 ) (2) -o (3) -io (2) o (2) -I (2) fo (2) -o (2) -a (2) -o (2) +o (2) 

i ymm,ny t t,n ytm,tm 2Htt ytt,mm <2,ymn I y mn u yml.ln ynl,lm ytt,mn ^ yil,mn 



(104) 
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The Einstein tensor components (|10[) are 



rr(2) 
"it 


= 4 2) 


- hR {2) - 


i (2) + 


i (2) 








- hR {A) - 


U 2) ^ (2) = 






r»(3) 


d(3) 
- W ti 


- ± (3) 


_i„( 2 ) _ 

2 •3im,mt 


- ±a (3) 4 


- ± (2) 


r»(2) 


d(2) 


+ %i?( 2 ) = 


= ±o (2) - 

2 &%3,mm 


i fl (2) .- 


V 2) . 



± (2) +io (2) 



By assuming the harmonic gauge 



(105) 



(2) _ (2) , (2) 

9tt,mm 9nn,mm ~r 9rnn,mn 



9 pa ^ = 

it is possible to simplify the components of Ricci tensor (|103[) . In fact for [i = one has 



(106) 



2fl <T T t sa n w - 2o w 4- rv w — 



,(2) _ ,„(3) j. (2) 

ymm. t 

and for fi = i 

z .y 1 err ~ 9tt,i T A 9 m i,m 9mm,i ~ u 

Differentiating (|10T[) with respect to t, a; J and (|108[) and with respect to t, one obtains 



(107) 



(108) 



(2) _ 2 (3) (2) 
J/tt.it ^9tm,mt t y mm tt ~ u 



(2) _ 2 (3) (2) 



(2) 



2o (2) - o (2) 

^rni.tm -Jrara.ti 



On the other side, combining (|110j) and (|111[) , we get 



(2) _ (2) (3) 
9mm,ti 9mi,tm 9mt,mi 



Finally, differentiating (|108p with respect to x 1 , one has : 



= 



(109) 
(110) 
(111) 

(112) 



o (2) +2 (2) 

9tt,ij ' z 9n 



9 



(2) 







and redefining indexes as j — > i, i —> j since these are mute indexes, we get 



(113) 



(2) (2) (2) 

ytt,iy 1 cfmy,im jmm.tj 



(114) 



The gauge transformation is h^ v = h^ v — f „ v ~ Cv,u when we perform a coordinate transformation as x'^ = + £ M with 0(£ 2 )<S 1. 
To obtain our gauge and the validity of the field equations for both perturbation and the vector £ M has to satisfy the harmonic 
condition = 0. 



23 



Combining (|Tl3t and (ITT4)) . we obtain 



.(2) , „(2) 



.,(2) 



(2) 
ymm,ij 



Relations (1109[) . (1112[) . (|115p guarantee to rewrite the Ricci tensor components (| 103[) at higher orders as 



n(2)| _ 1a (2 



(2) 



(4) 



1 (2) (2) _ 1 (2) 
l9rang t t.mn 29tt,tt 



Rf^HG 



(2) 
2 "fij' 



and the Ricci scalar (| 104)) becomes 



i?( 2 )| HG = |A^f-iA^ 



(115) 



(116) 



(117) 



where V and A are, respectively, the gradient and the Laplacian in flat space. The Einstein tensor components (|10[) 
in the harmonic gauge are 



( G[f\ HG = lA 9 ^ + lA 9 ,. 



Wt l#G 



^(3) I 



^(2) I 
Cry I^G 



(2) , 1A„(2) 



(3) 



1 Ao (2) + 

2 aij ~ 4 



(118) 



Ao (2) - Ao (2) 



On the matter side, i.e. right-hand side of the field equations ([9]), we start with the general definition of the 
energy-momentum tensor of a perfect fluid ([7]) with additional energy II 



Tfj, u = (p + Up + p)u^u v -pg^ 
The explicit form of the energy-momentum tensor can be derived as follows 



(119) 



T tt = P + p{v 2 -2U + Tl)+p 



v 2 \ £+v 2 + 2V + 11) n 2.17 P 



-pv' + p 



-r-\ ^ + 2v + v 2 + n )+h t 



i V 



Tij = pv l v j + p 



»V n + 2 + 4F + u 2 + 2/7 - 2« c 5 c(i %) + 2^% 



(120) 



As a first application of these results, let us now take into account the simplest case, that is GR. 
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A. The Newtonian and Post-Newtonian Limit of General Relativity 

Einstein Eqs. ([9]) can be rewritten as 



Ru,v — X 



T 



(121) 



From the interpretation of stress-energy tensor components as energy density, momentum density and momentum 
flux, we have T tt , T t i and at the various orders 



T tt = 4 0) +T^+O(4) 



-.(2) 



T tt = T t { p + 0(3) 



(122) 



rp . rp(2) 

J ij - 1 ij 



0(4) 



where denotes the term in T M „ of order M/r 3 v N . In particular T^"' is the density of rest- mass, while T^' is 

the non-relativistic part of the energy density. What we need is the tensor 



,(0) 



.(2) 



T 



(123) 



GM/f is of order v 2 , so O and (TT2"2")) give 



s ti 

Sj-j = 



q(0) , o(2) 
°tt "+" D tt 



0(6) 



S«+0(3) 



(124) 



4V 0(2) 



where <S^^ denotes the term in S^,, of order M/r 3 v N . In particular 



c(0) _ 1 T (0) 
J tt — 2 tt 



d(2) 



l7i(2) , 171(2) 
2 1 U T 2 Jmm 



(i) 



(125) 



o(°) _ I5..7" 



(0) 



Substituting Eqs. (|116[) and (|124[) in Eqs. (|121j) . we find that the field equations in harmonic coordinates are indeed 
consistent with the expansions we are using, and give 



(2) 



R 



p(4) 
H tt 



R 



(3) 



(0) 



XS, 



xsg> 



(0) 



(126) 



R 



(2) 



XS, 



(0) 
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and, in particular, 



A 5t ( t 2) = 



A 



-,(2) 



Tt { t ] + T, 



(2) 
mm 



(2) (2) (2) 
9mn9tt,mn "+" 5tt,M "+" 



v<4 2) l 2 



Agf = 2A-T f 



(i) 



(0) 



(127) 



From the first one of (|127[) . we find, as expected, the Newtonian result: 



(2) 

9tt 



-— / d J x 

477/ |X — X' 



2G / ef x' 



2$(x) 



(128) 



where ^(x) is the gravitational potential. In fact if we consider a point-like source (|5Tj) we find 



$(x) = - 



GAf 



(129) 



From the third and fourth equations of system (|127|) . we find that 



9t 



(3) 



= -2l/^x'^?^(x) 



The second equation of (|127p can be rewritten as follows 



(130) 



A 



(4) 

9tt 



2$' 



= X 



T (2) , T (2) 
t£ ^ mm 



- 8$A$ + 2$ , 



(131) 



and the solution for gffl is 



g $> = 2 ** - £ r^T^ni^y) 2 r *wa«*v) j_ r^^o = 2T(x) (132) 

y " 47T / X-X' 7T / x-x' 2?r / x-x' v ; v y 



.(2) 



By using the equations at second order, we obtain the final expression for the correction at fourth order in the 
time-time component of the metric: 



T(x) = ^-£/^^^ + ^ + f ^ x ,*(x')T£V) 1 fl2 *(*) 



4tt x-x' 



(133) 



We can rewrite the metric expression (|98|) as follows 



1 + 2$ + 2 T Z T 

Z _^.(l_ 2 $) y 



(134) 



where Z are higher order terms that can be assumed null at this approximation level. 
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Finally the Lagrangian of a particle in presence of a gravitational field can be expressed as proportional to the 
invariant distance ds 1 / 2 , thus we have: 



L = \g P a 



dx p dx° 



a\ 1/2 



dt dt 



= [gtt + 2gt m v r ' 



g mn v m v n ^j 1 = (l + gW + gW+2gWv m -v 2 +g£lv m v^ ' (135) 



which, to the 0(2) order, reduces to the classic Newtonian Lagrangian of a test particle L 



New 



= 1 + 2$ - 



1/2 



wher 



dt 



and |v| 2 = v m v m . As matter of fact, post-Newtonian physics has to involve higher than 0(2) order 



terms in the Lagrangian. In fact we obtain 



1 



T 



(136) 



An important remark concerns the odd-order perturbation terms 0(1) or 0(3). Since, these terms contain odd 
powers of velocity v or of time derivatives, they are related to the energy dissipation or absorption by the system. 
Nevertheless, the mass-energy conservation prevents the energy and mass losses and, as a consequence, prevents, in the 
Newtonian limit, terms of 0(1) and 0(3) orders in the Lagrangian. If one takes into account contributions higher than 
0(4) order, different theories give different predictions. GR, for example, due to the conservation of post-Newtonian 
energy, forbids terms of 0(5) order; on the other hand, terms of 0(7) order can appear and are related to the energy 
lost by means of the gravitational radiation. 



V. THE NEWTONIAN LIMIT OF /(7?)-GRAVITY BY THE O'HANLON THEROY ANALOGY 

Let us start our analysis of Newtonian and post-Newtonian limits of Extended Theories of Gravity discussing the 
possible shortcomings related to the use of analogies in the weak field limit approximation. As briefly pointed out 
above, some authors have claimed that FOG models are characterized by an ill defined behavior in the Newtonian 
regime. In particular, in a series of papers [l?} it is addressed that post-Newtonian corrections of the gravitational 
potential violate experimental constraints since these quantities can be recovered by a direct analogy with Brans- 
Dicke Gravity 50] simply supposing the Brans-Dicke parameter ujbd in Eq. (I2~5j) vanishing for /(i?)-gravity. Actually 
, despite the calculations of the Newtonian and the post-Newtonian limit of /(i?)-gravity, performed in a rigorous 
manner, have showed that this is not the case [44l I51M531 ]. it remains to clarify why the analogy with Brans-Dicke 
Gravity seems to fail its predictions. The issue is easily overcame once the correct analogy between /(i?)-gravity and 
the Brans-Dicke theory is taken into account. 

It can be easily shown that, /(-R)-gravity models can be rewritten in term of a scalar-field Lagrangian non minimally 
coupled with gravity but without any kinetic term implying ujbd = 0. Actually, the simplest case of Scalar- Tensor 
Gravity models has been introduced some decades ago by Brans and Dicke in order to give a general mechanism 
capable of explaining the inertial forces by means of a background gravitational interaction. The explicit expression 
of such gravitational action is (|25l) . while the general action of /(-R)-gravity is (fT5|) when f(X,Y,Z,) = f(R). As 
said above, /(i?)-gravity can be recast as a Scalar- Tensor theory by introducing a suitable scalar field <fi which non- 
minimally couples with the gravity sector. It is important to remark that such an analogy holds in a formalism in 
which the scalar field displays no kinetic term but it is characterized by means of a self-interaction potential which 
determines the whole dynamics (O'H anion Lagrangian) (Bij . We can resume the actions as follow 



A f P=fd*xV=9 



f(R) + XC r , 



cpR-io BD ^-+X£ ri 



(137) 



A°» = fd*xV=9 



4>R + V{4>) + XC r , 



where JF means that we are considering all theories in the Jordan frame. This consideration, therefore, implies that 
the scalar field Lagrangian equivalent to the purely geometrical /(i?)-gravity turns out to be quite different with 
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respect to the ordinary Brans-Dicke action in (|25j) . This point represents a crucial aspect of our analysis. In fact, 
as we will show, such a difference implies completely different results in the Newtonian limit of the two models and, 
consequently, shows that it is misleading to extend predictions from the PPN approximation of Brans-Dicke models 
to /(i?)-gravity. Considering natural units, the O'Hanlon Lagrangian [54| is the third of (|137l) . The field equations 
are obtained by varying the action with respect to both g^ v and <f> which now represent the dynamical variables (the 
same field equations are given setting uj(<fi) = and F{<f>) — <fi in Eqs. ([26)) ). Thus, one obtains 



' fyG^v - ^Vi^g^ - <f> ;)JV +g llv U(j) = XT^ U 

R , dV{4>) 







(138) 



< deb 3 



v{<t>) 



X 



T 



where the second line of (|1 38[) is the field equation for cf>. While the third equation is a combination of the trace of 
the first and the second ones. Field equations for /(i?)-gravity are obtained from (fT8|) 



< 

^ 3D/' + f'R -2/ = XT 
The two approaches can be mapped one into the other considering the following equivalences 



(139) 



'</> = /' 

< V(<(>) = f-f'R 



(140) 



;dV{<t>) 



2V{4>) = f'R-2f 



where the Jacobian matrix of the transformation <p ■<=>■ /' has to be non- vanishing. Henceforth we can consider 
instead of (|139[) a new set of field equations determined by the equivalence of /(iZ)-gravity with the O'Hanlon approach 

M 



<t>R^ + ± [ V{4>) + )g» v - = XT, 



(141) 



3 



2V{cj>) 



x 



T 



where T^ = T^„ — ^Tg^. Let us, now, calculate the Newtonian limit of field equations (|14ip . We take into 
account the perturbations of metric tensor g^ v in Eqs. (l98|) up to C(2)-order and also for scalar field an analogous 
perturbation with respect to the background value 



+ ■ 



(142) 



The differential operators turn out to be approximated as 



□ « df - A and V^V^ d 



f.LV 



(143) 



Actually in order to simplify calculations we can exploit the intrinsic gauge freedom in the metric definition. In 
particular, we choose the harmonic gauge (|106l) and the expressions of Ricci tensor components are given by (|116l) . 
In relation to the adopted approximation we coherently develop the self-interaction potential at second order. In 
particular, the quantities in (|141[) read 
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^VM _ 2F( 0) * _ 2 V(0(°))+ 



d4> 2 d<j> 



0( 2 ) 



The field equations (|14ip . solved at 0(O)-order of approximation, provide the two solutions 



(144) 



Vm=0 and ^1=0 



(145) 



which fix the O(0)-order terms in the development of the self-interaction potential; therefore we have 



V{4>) + 0^ * ^) d ^fl^) = 3m 2 0( 2 ) 
0^M - 2V(<P) ~ ^^1^) = 3m 2 0( 2 ) 



(146) 



where constant factors 0(°) d — - have been condensed within the quantity 3m 27 . Such a constant can be easily 
interpreted as a mass term as will become clearer in the following. Now, taking into account the above simplifications, 
we can rewrite field Eqs. (|141[) at the at C(2)-order in the form 
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The scalar field solution can be easily obtained from Eq. (|149l) as 



(148) 
(149) 



X r d 3 k f(°)(k)e lk - 5 
3" J (2tt) 3 / 2 k 2 + to 2 

(2) 



where T(°)(k) is the Fourier transform of the trace T(°). While for gffl and <jrj?' we have 



(150) 
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m 2 /> x ^ (2) (x') 



4tt0(o) 



(151) 
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y (0) fx ,N i 
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4tt0(o) 
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|x - X' 
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(152) 



7 The factor 3 is introduced to simplify an analogous factor present in the field equations 11411 . 
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The above three solutions represent a completely general result. In particular adopting transformations (I140p . one 
can straightforwardly obtain the solutions in the /(-R)-scheme. 

Let us analyze the above results with an example. We can consider a FOG Lagrangian of the form 



f(R) = ai R + a 2 R 2 (153) 

so that the scalar field reads <f> — a± + 2a 2 R (ai and a 2 are arbitrary constants). The self-interaction potential turns 

out the be V((j>) — — satisfying the conditions V{a\) = and V'(ai) = 0. In relation with the definition of 

the scalar field, we can opportunely identify a\ with a constant value (f>( ' — a\. Furthermore, the scalar field "mass" 
can be expressed in term of the Lagrangian parameters as follows 



Since the Ricci scalar at the second order reads 



- 2 = 3^ 0) T = "S (154) 



n ~ = = A / rf 3 k f(°)(k)e' k x 
2a 2 6a 2 i (2tt) 3 / 2 k 2 + m 2 

if we consider a point-like source (|51[) therefore we obtain 



(155) 



GM /• 3 e lk x _ nfr g m 2 e" m l x l 
k 2 + m 2 ~ "V 2~ ji[ 

The immediate consequence is that the solution for the scalar field <f> at second order is 

nr a e-" l l x l 



= 2 -« (2 > - Vif V (157) 

while the complete scalar field solution up to the second order of perturbation is given by 

/— p -m\x\ 

rf^r (158) 

Once the behavior of the scalar field has been obtained up to the second order of perturbation, in the same way, one 

can deduce the expressions for g$ and gyj\ where T,^' = |p and E^- 1 = \p5ij — £y. As matter of fact the 
metric solutions at the second order of perturbation are 
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This quantity, which is directly related to the gravitational potential, shows that the gravitational potential of the 
O'Hanlon Lagrangian is non-Newtonian. Such a behavior prevents from obtaining a natural definition of the PPN 
parameters as corrections to the Newtonian potential. As matter of fact, since it is indeed not true that a generic f{R)- 
gravity model corresponds to a Brans-Dicke model with ojbd = coherently to its post- Newtonian approximation. 
In particular it turns out to be wrong considering the PPN parameter 7 = 2+Zbd ( see ' ^ or exam P^ e i UM) °f Brans - 
Dicke gravity and evaluating this at ojbd = so that one gets 7 = 1/2 as derived in [47$ . 
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Differently, because of the presence of the self-interaction potential V(4>), in the O'Hanlon Lagrangian, a Yukawa 
like correction appears in the Newtonian limit. As matter of fact, one obtains a different gravitational potential with 
respect to the standard Newtonian one and, as matter of fact, the fourth order corrections in term of the v/c ratio 
(Newtonian level), have to be evaluated in a complete new way. In other words, considering a Brans-Dicke Lagrangian 
and an O'Hanlon one, despite their similar structure, implies different predictions in the weak field and small velocity 
limits. Such a result represents a significant argument against the claim that FOG models can be ruled out only on 
the bases of the analogy with Brans-Dicke PPN parameters. 

An important consideration is in order at this point on the meaning of PPN-parameters 7 and /3, defined as a 
correction to the Newtonian-like behavior of the gravitational potentials (|69|) . Actually, if we consider the limit 
f(R) — > R, from (I159[) and set a± = 2/3 {a\ is an arbitrary constant), we have 
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which suggest that the PPN parameter 7, in this limit, results 1/2 which is in striking contrast with GR predictions 
(7 ~ 1). Such a result is not surprising. In fact, the GR limit of the O'Hanlon Lagrangian requires <f> ~ const and 
V{4>) —¥ but such approximations induce mathematical inconsistencies in the field equations of f(R) -gravity, once 
these have been obtained by a given O'Hanlon Lagrangian. Actually this is a general issue of O'Hanlon Lagrangian. 
In fact it can be demonstrated that field Eqs. (|141l) do not reduce to the standard GR ones (for V{<f) — > and 
<t> const) since we have 



</>iV + i (v(<t>) + 0^ V„ - ^ v = xz^ 

In fact £ M „ components read E tt = |p and = \p&ij — §Ett $ij in place of Stt = \p and Sij = \p^ij = Stt $ij 
as it should be for the GR field Eqs. (1121[) . Such a pathology emerges also when the GR limit is performed from a 
pure Brans-Dicke Lagrangian. In such a case, in order to match the Hilbert-Einstein Lagrangian, one needs <j> ~ const 
and ojbd = 0, the immediate consequence is that the PPN parameter 7 turns out to be 1/2, while it is well known 
that Brans-Dicke model fulfils low energy limit prescriptions in the limit U) — > 00. Even in this case, the problem, 
with respect to the GR predictions, is that the GR limit of the model introduces inconsistencies in the field equations. 
In other words, it is not possible to impose the same transformation which leads the Brans-Dicke theory into GR at 
the Lagrangian level on the solutions obtained by solving the field equations descending from the general Lagrangian. 
The relevant aspect of this discussion is that considering a /(i?)-model, in analogy with the O'Hanlon Lagrangian 
and supposing that the self-interaction potential is negligible, introduces a pathological behaviour in dynamics which 
results in obtaining a PPN parameter 7 = 1/2. This is what happens when an effective approximation scheme is 
introduced in the field equations in order to calculate the weak field limit of FOG by means of Brans-Dicke model. 
Such a result seems, from another point of view, to enforce the claim that FOG models have to be carefully investigated 
in this limit and their analogy with scalar-tensor gravity should be opportunely considered. 
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A. Scalar- Tensor Gravity in Jordan and Einstein frames 



Up to now, we have discussed the weak field and small velocity limit of FOG in term of Brans-Dicke like La- 
grangian remaining in the Jordan frame. We have considered the weak field and small velocity limit when a conformal 
transformation (f27|) is applied to the O'Hanlon Lagrangian. Now we want to analyze the differences of considering 
a Scalar- Tensor theory in the Jordan frame and in the Einstein frame. The Scalar- Tensor action Ajp in the Jordan 
frame (j2~i)) is linked to the action A^f F in the Einstein frame (|3T)]) via the transformations (|3"Tj) between the quantities 
in the two frames. The O'Hanlon theory in the Jordan frame is recovered from Eq. (|24p imposing F(cf>) — <f> and 
uj(4>) = 0. Action (j30]), in the Einstein frame results simplified and the transformation between the two scalar fields 
reads 
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2 3A# 
S2(^)d^ = — ^--7o- 
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If we suppose £l(p) = — f2o < we have 
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where A = \ I and /c is an integration constant. The O'Hanlon theory, transformed in the Einstein frame, is 
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The field equations are 



A 2 vl2 /A 



(164) 



1 Al e T2A V1 



2O fi^+^eT 2 V[^( fce ±A V)T 2Ay(fc e ± ^)]+A'£ m ^ = 



7? — — 2Lf"P 4- ^-lo w' c 
{ Jl — 2A T A r^af 



2A, 



T 2A * , y(fe ±A «') 



(165) 



where the matter tensor, which now coupled with the scalar field ip, in the Einstein frame reads 
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and 
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Actually, in order to calculate the weak field and small velocity limit of the model in the Einstein frame, we can 



develop the two scalar fields at the second order 
value. This choice gives the relations 



( 2 ) and ip ~ <p(°> + <p( 2 ) with respect to a background 
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Let us consider the conformal transformation g pi/ — j^g^v From this relation, considering Eq. (|163l) . one obtains for 
= A 
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As matter of fact, since g^ ] = 2$ Ji? , g^f = 2* JF % and g£ ] = 2<5> EF , g^ ] = 2 V EF <% from Eqs-dUlD, relevant 
relations emerge linking the gravitational potentials between Jordan and Einstein frames 
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20(°) ~ * + 2^ 

If we introduce the variations of two potentials, A4> = $ ,7F — $> EF and ASP = ^> ,IF — ^> EF , we obtain 



A$ = -A* = = <x a 2 cx /"(0) 



(171) 



here specified in the case of Lagrangian (|153l) . 

Prom the above expressions, one can notice that there is an evident difference between the behavior of the two 
gravitational potentials in the two frames |55l |. Such a result suggests that, at the Newtonian level, it is possible 
to discriminate between the two frames thus one can deduce what is the true physical one. In particular, once, the 
gravitational potential is calculated in the Jordan frame and the dynamical evolution of 4> is taken into account at the 
suitable perturbation level, this can be substituted in the first of (|170l) to obtain its Einstein frame evolution. The 
final step is that the two potentials have to be matched with experimental data in order to select what is the true 
physical solution. 



VI. THE NEWTONIAN LIMIT OF /(^-GRAVITY IN STANDARD COORDINATES 

The Newtonian limit of FOG can be worked out by comparing its viability with respect to the standard results 
of GR. Here, we investigate the limit in the metric approach, refraining from exploiting the formal equivalence of 
FOG with specific Scalar- Tensor theories, i.e. we work in the Jordan frame in order to avoid possible misleading 
interpretations of the results [45[ . 

Considering the Taylor expansion of a generic /(-R)-gravity model, it is possible to obtain general solutions in term 
of the metric coefficients up to the third order of approximation. Furthermore, it is possible to show that the Birkhoff 
theorem is not a general result for /(i?)-gravity since time-dependent evolution of spherically symmetric solutions can 
be achieved depending on the order of perturbations. 

Exploiting the formalism of Newtonian and post-Newtonian approximations previously described, we can develop 
a systematic analysis in the limits of weak field and small velocities for /(i?)-gravity. We are going to assume, as 
background, a spherically symmetric spacetime and we are going to investigate the vacuum case. Considering the 
metric (jl2")l , we ha,ve, for a. given 
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The determinant reads 
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Christoffel symbols (|102[) are 
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while the Ricci tensor components (| 103[) are 
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and, finally, the Ricci scalar expression is 
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By metric tensor (|172p and by inserting it into the field equations (I139p . one obtains 



(177) 



Hfjj/ / R^v 2 /-9t" y / I ^uvRj I" 'uvR,r g^v 



g' 



g rr In ^ r i?,r + + <T"-K,r 



R,ii.R,u — g^v [ g u R,t 2 + g 7r R,r 2 



H = f'R - 2/ + 3/" 
+3/"' 



g u R/+g rr R/ 



</» + <?" In ^ J iJ, t + <T r + <T In ^=5. fl,r + ff"iJ,tt + g rr R, 



(178) 



In order to derive the Newtonian and post-Newtonian approximations for a generic /(i?)-function , one should 
specify the f(R) - Lagrangian into the field Eqs. (|178p . This is a crucial point because once a certain Lagrangian 
is chosen, one will obtain a particular approximation referred to such a choice. This means to lose any general 
prescription and to obtain corrections to the Newtonian potential, $(x), which refer "univocally" to the considered 
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/(i?)-function. Alternatively, one can restrict to analytic /(i?)-functions expandable with respect to a certain value 
R = R = constant or, at least, its non-analytic part, if exists at all, goes to zero faster than R n , with n > 2 
at R — >• 0. In general, such theories are physically interesting and allow to recover the GR results and the correct 
boundary and asymptotic conditions. Then we assume 



f(R) = E LA j 2l (R - Ro) n ^fo + hR + hR 2 ■ f:jr • ... (179) 

n 

One has to note that the expansion (|179p . also if similar to ([75]) . presents some differences. In fact R^ is a general 
space-time function linked to the background metric gffi given in (|74[) . Here i?o is a constant value of scalar curvature, 
which can be negligible in weak field limit approximation. Besides, the coefficients fo, /i, fi, /3 are not proportional, 
respectively, to zero-th, first, second, third coefficient of the Taylor expansion of f(R). In fact, we have 

' fo = f(Ro) - Rof'(Ro) + k R of"(Ro) - lR 3 of"'(Ro) 
h = f'(Ro) - Rof'(Ro) + kRof"'(Ro) 

(180) 

fi = ^f'(Ro) - \Rof'"{Ro) 

k h = y"'(Ro) 

If we consider a flat background, then Ro = and the coefficients fo, /1, fa, are the terms of Taylor series. 
But if we are searching for solutions at Newtonian and (possibility) post-Newtonian level, we have to consider a 
vanishing background scalar curvature. It is possible to obtain the Newtonian and post-Newtonian approximation of 
/(ii)-gravity considering such an expansion (|179l) into the field Eqs. (| 1 78[) and to expand the system up to the orders 
O(0), 0(2), 0(3) and 0(4). This approach provides general results and specific (analytic) Lagrangians are selected 
by the coefficients fi in Eq. (|179l) . Developing the equations in the case of vanishing matter, i.e. T^ v = 0, we have 
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and, in particular, from the O(0) order approximation, one obtains 



fo = (182) 

which trivially follows from the above assumption that the space-time is asymptotically Minkowski (asymptotically 
flat background) . This result suggests a first consideration. If the Lagrangian is expanded around a vanishing value 
of the Ricci scalar (R — 0), the relation implies that the cosmological constant contribution has to be zero 

whatever is the f{R)- gravity model. 
If we now consider the C(2)-order approximation, system (|181[) . in the vacuum case, results to be 
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The last equation of the system l|183|) is the definition of Ricci scalar ((5|) at 0(2)-order. The trace equation (the 
fourth line in Eqs. (|183|) ). in particular, provides a differential equation with respect to the Ricci scalar which allows 
to solve, if sign[/i] = — sign[/ 2 ], the system (|183[) at C(2)-order. The solutions are 
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with the dimension of length" 1 . We note that the definition of mass (|185[) is compatible with the definition of (|154|) . 
Let us notice that the integration constant So has to be dimensionless, Si has the dimension of length, while the 
time - dependent functions 82 and £3, respectively, have the dimensions of length -1 and length" 2 . The functions 
8i(t) (i = 2, 3) are completely arbitrary since the differential equation system (|183p contains only spatial derivatives. 
Besides, the integration constant So can be set to zero, as in the theory of the potential, since it represents an 
unessential additive quantity. When wc consider the limit f(R) — > R, in the case of a point-like source (f5Tj) , we 
recover the perturbed version of standard Schwarzschild solution (|5^|) at C(2)-order with <5i = r g . In order to match 
at infinity the Minkowskian prescription of the metric, we discard the Yukawa growing mode present in (|184[) . Then 
we have, in standard coordinates, 
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At this point one can provide the solution in term of the gravitational potential. In such a case, we have an exp licit 
Newtonian- like term into the definition, according to previous results obtained with less rigorous methods [l2l l56l|. 
The first of (I184p provides the second order solution in term of the metric expansion (see the definition (|172p ). but, 
this term coincides with the gravitational potential at the Newtonian order (|128[) . In particular the gravitational 
potential of a FOG-model, analytic in the Ricci scalar R, is 
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As first remark, one has to notice that the structure of the potential (|187|) . for a given /(i?)-gravity, is determined 
only by the parameter m, (|185p . which depends on the first and the second derivative of the /(i?)-function, once 
developed around a vanishing value of the Ricci scalar. 
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As second remark, the solution for the gravitational potential has a Yukawa-like behavior depending on a charac- 
teristic length on which it evolves. 

In other words, the correction to the Newtonian gravitational potential is always characterized by a Yukawa-like 
correction and only the first two terms of the Taylor expansion of a generic f (i?) -function turn out to be relevant. 
This is indeed a general result. 

Let us now consider system (|181j) at third order approximation. The first important issue is that, at this order, one 
has to consider even the off-diagonal equation 



fig% + 2f 2 rR$ = (f 88) 

(2) 

which relates the time derivative of the Ricci scalar to the time derivative of gW ■ From this relation, it is possible 
to draw a relevant conclusion. One can deduce that, if the Ricci scalar depends on time, so it is the same for the 
metric components and even the gravitational potential turns out to be time-dependent. This result agrees with the 
analysis provided in (3(| where a complete description of the weak field limit of FOG has been provided in term of 
the dynamical evolution of the Ricci scalar. Moreover it has been demonstrated that supposing the time independent 
Ricci scalar, a static spherically symmetric solution is found. 

Eq. (jf 88[) confirms this result and provides the formal theoretical explanation of such a behavior. In particular, 
together with Eqs. (|I86[) . it suggests that if one considers the problem at a lower level of approximation (i.e. the 
second order) the background spacetime metric can be only factorized with a function space-depending and an 
arbitrary function time-depending. Then the Birkhoff theorem at Newtonian level is modified. The static solutions 
according to the Birkhoff theorem in GR are not directly obtained. Obviously this is still no more verified when 
the problem is faced with approximations of higher order. In other words, the debated issue to prove the validity of 
the Birkhoff theorem in FOG, finds here its physical answer. In (36|, the validity of this theorem is demonstrated 
for FOG only when the Ricci scalar is time-independent or, in addition, when the solutions of (II 391) are investigated 
up to the second order of approximation in the metric coefficients (|I72|) . Therefore, the Birkhoff theorem is not a 
general result for FOG but, on the other hand, in the limit of small velocities and weak fields (which is enough to deal 
with the Solar System gravitational experiments), one can assume that the gravitational potential is effectively time 
independent according to U86\) and {18ty . 

The above results fix a fundamental difference between GR and FOG. While, in GR, a spherically symmetric solution 
represents a stationary and static configuration difficult to be related to a cosmological background evolution, this 
is no more true in the case of generic FOGs. In the latter case, a spherically symmetric background can show time- 
dependent evolution together with the radial dependence. In this sense, a relation between a spherical solution and 
the cosmological Hubble flow could be, in principle, achieved. 



VII. THE NEWTONIAN AND POST-NEWTONIAN LIMIT OF /(i?)-GRAVITY IN ISOTROPIC 

COORDINATES 

Let us consider now a form of metric tensor generalizing the metric (|I72[) . It is interesting, using the isotropic 
coordinates (|4I[) and a more general approach, to solve the field equations as shown in 46]. The metric which we 
take into account is the following 

(l+g£ ) (t, X )+git\t,x) + ... g£\t, X ) + ... \ 
^ [ + -^ +fl g>(t,x) + ...J (189) 

and the set of coordinates adopted is a;'' = (t,x) = (t, x 1 , x 2 , x 3 ). The n-th derivative of function / can be developed 
as in Eq. (fT?9")) with R = 

f n (R) ~ f n {R {2) + i? (4) + . . . ) ~ /™(0) + f n+1 (0)R {2) + / Il+1 (0)i? (4) + i/" +2 (0)i? (2)2 + . . . (190) 

From lowest order of field Eqs. (|72")l , we find the same condition (|I82[) (f(0) — 0). At 0(2)-order (Newtonian level) 
Eqs. (|72|) becomes 
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t ff(2) = -3/"(0)Ai?( 2 ) - f(0)R^ = XTM 
while at 0(3)-order becomes 
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Remembering the expressions of Christoffel symbols and using the following approximation for the determinant of 
metric tensor In y/— g ~ ^[9tt ~ 9mm] + • • ■ , at 0(4)-order, we have 
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Note that the propagation of Ricci scalar i?( 4 ) has the same dynamics of previous one (third line of Eqs. (|19ip ). The 
complete knowledge of correction at fourth order for the tt-component of Ricci tensor fix the third derivative of f(R) 
in R = 0. Also at this level, there is a degeneracy of /(iZ)-theory: different theories, considering only the first three 
derivatives, admit the same gravitational field without radiation emission. 

We want to rewrite and generalize the outcome of Eqs. (|186|) by introducing the Green function method (we 
remember that the Newtonian limit corresponds also to the linearization of field equations). Let us start from the 
trace equation. The solution for the Ricci scalar R^ in the third line of Eqs. (|191l) is 
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where we defined 
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and <?(x, x') is the Green function of field operator A — m 2 . We note that the definition of mass (|195[) is compatible 
with the definitions of (I154[) and (|185l) . 



The solution for g\*\ from the first line of (|191l) by considering that R^ = hAg^ (first line of (I103[) or (jllGI) 1 ). is 
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We can check immediately that when f{R) — >• R we find g^ (t, x) — > — 2G J d 3 x' j^^j ■ The expression (|196[) is 

the "modified" gravitational potential (here we have a factor 2) for /(i?)-gravity and generalize, for any matter 
distribution, the outcome (|187[) . This solution, which is the Newtonian limit of f{R)- gravity, is also gauge-free. 

Since we have a linearized version of field equations, this limit corresponds to one of the Einstein equation and the 
linear superposition is satisfied. So the tt-component of energy-momentum tensor is, in this limit, the sum of mass- 



energy- volume density of sources, that is: T t 
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£ a M a <5(x — x a ) where S(x) is the delta function. 



As it is evident, the Gauss theorem is not valid since the force law is not cx |x| . The equivalence between a 
spherically symmetric distribution and a point-like distribution is not valid and how the matter is distributed in the 
space becomes extremely important in this situation. However, we have to say that the Bianchi identities hold in any 
case so the consistency of the theory is guaranteed. 

From field Eq. (|192[) . by using the gauge harmonic condition (|106[) . we find the general solution for g!jjp 
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_______ ^ ^ \ 

The choice of harmonic gauge enable us to solve Eq. (|192j) but we lose information on time evolution of g tt (f, x). 
This is important to obtain, at least in perturbative approach, some information about the Birkhoff theorem. By 
hypothesizing a perturbative approach (Newtonian- like), we confine time-evolution only to the timel- variation of 
matter source. In fact in this working hypothesis, the motion of bodies embedded in gravitational fields evolves very 
slowly with respect to the internal motions of matter. Then we have, in any case, an instantaneous readjustment of 
spacetime. In other words, the motion of bodies is adiabatic and it enables us to factorize the solutions that, by a 
time transformation, become static solutions. 

Still more, also the corrections to the gravito-magnetic effects (|192l) are depending on the only first two derivatives 
of f(R) in R = 0. This means that different theories, from the third derivative on, admit the same Newtonian 
solution. 

_____ ______ f<2\ 

From second line of (|191[) . by using the gauge harmonic condition (jl06[) . the solution for g\-' follows 
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where J7| x | represents the integration volume with radius |x| (for the details see [26|). By the solutions (|196l) . (1197[) . 
()198[1 we can affirm that it is possible to have solution non-Ricci-flat in vacuum. This means that: Higher Order 
Gravity mimics a matter source. It is evident, from (|196|) . that the Ricci scalar can be considered a "matter source" 
which curves the spacetime also in absence of ordinary matter. Then it is clear also that the knowledge of behavior 
of Ricci scalar inside mass distribution is fundamental to obtain the behavior of metric tensor outside the matter. 

From the fourth order of field equations, we note also that the Ricci scalar R^ propagates with the same m (the 
second line of (|193p ) and the solution at second order originates a supplementary matter source in r.h.s. of (f72"|) . The 
solution is 
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Also in this case we can have a non- vanishing curvature in absence of matter. The solution for g[f \ from the first 
line of (M]), is 
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(t,x')g$ mn (ty) ~ g%{t,x>) - |V x < 5( (2 V,x')| 2 - R^(t, X ') - 5t (2) (t,x')i? (2) (*,x') 

i?( 2 ) 2 (i,x') J?( 2 )(t,xOA^ 2) (t,x') 



67rm^ 



+ £L m (*> x')i? (2) ,»(*, x') + AflW (t, x') 



+ 5t (2) (t,x')A J R( 2 )(t,x') +^)(t,x')i? (2) , m „(t,x') - lv 5 ( n l(t,x') • Vi?( 2 )(t,x') 



(201) 



In summary, we have shown the more general solutions of field equations of /(i?)-gravity in the Newtonian and 
post-Newtonian limits assuming a coordinates transformation where the gauge harmonic condition is verified. Now 
we shall apply such an approach to obtain the explicit form of the metric tensor for a static and spherically symmetric 
matter source. 

A. Solutions generated by an extended spherically symmetric source with harmonic gauge conditions 



Let us consider a spherical source with mass M and radius £. Since the metric is given by Eqs. ()189p . the energy- 
momentum tensor (|5ip becomes 



T«(t,x) ~ p(pc) + p{x)g%\t,x) = T t ( t 0) (t,x)+T t (2) (t,x) 



t = P (x) = r(°>(t,x) 



(202) 



where the density is 



^ (x) = |x|) 



(203) 



O(x) is the Heaviside function. We are not interested here to the internal structure of the source. The possible choices 
of the Green functions of the field operator A — to 2 , for spherically symmetric systems {i.e. C/(x, x') = <?(|x — x'|)), 
are the following 



1 e" 



<?(x,x') 



Ci 



47T |x — x'| 



if to 2 > 



(204) 



|x-x'| 



• 2 ^r if m 2 < 



with Ci + C 2 = 



In the Newtonian limit of GR, the equation for the gravitational potential, generated by a point-like source 



AxSjVe-iD.mecft.Cx-x') = -47T<5(x-x') (205) 

is not satisfied by the Green functions (|204l) . If we consider the flux of gravitational field gNew.mech. defined as 



fsNew.mech. 



GM(x - x') 



,/|3 



— — GMV ^QNew.mech.^ ~ x ') 



(206) 



we obtain, as standard, the Gauss theorem 
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d£ gNew.mech. ■ h <X M 



(207) 



where £ is a generic two-dimensional surface and n its surface normal vector. The flux of field gNew.mech. ° n 

the 

surface S is proportional to the matter content M, inside the surface independently of the particular shape of surface 
(Gauss theorem, or Newton theorem for the gravitational field [13 )• On the other hand, if we consider the flux defined 
by the new Green function, its value is not proportional to the enclosed mass but depends on the particular choice of 
the surface 



rfS gNew.mech. ' ft OC M S 



(208) 



Hence Ms is a mass- function depending on the surface S. Then we have to find the solution inside/outside the matter 
distribution by evaluating the integral quantities and imposing the boundary condition on the separation surface. 
We have to note that, for any function of modulus ^(|x|), it is 



J d 3 x'e(x,x>(x') = -hj d \A\A 2 H\A) JJ d^' 



sin /9' p - m \/l x l 2 + l x 'l 2 ~ 2 l x ll x 'l COSQ 

rffl / S 6 (209) 
/,, v /|x| 2 + |x'| 2 -2|x||x'|cosa 



where cos a = cos 9 cos 9' + sin 6* sin 0' cos(</> — cf>') and a is the angle between two vectors x, x'. In the spherically 
symmetric case we can choose 9 = without losing generality (the symmetry of system is independent by the angle). 
By making the angular integration we get 



I = — 



1 



2to|x| 



d|x'||x'|/i(|x'|) 



e -m||x|-|x'|| _ e -m(|x| + |x'|) 



(210) 



An analogous relation is useful also for the Green function of Newtonian mechanics |x — x 



'i-i 



2tt 



J d|x'||x'| 



M|x'|) 



(211) 



1. Solutions at 0(2)- and (D(3)-order 
By supposing that m 2 > (i.e. sign[/'(0)] = — sign[/"(0)] (an analogous condition used in (I184p ). the Ricci scalar 

GEO; is 8 



i?< 2 '(i,x) = 



3r ( 



l-r m5 (l+mO 



sinh mlxl 



-m|x| 

- |x|) - r 9 rr?F{tO—— 6(|x| - 



(212) 



where we introduced a shape function 



F(x) = 3- 



i cosh mx — sinh ? 



(213) 



The solutions of (fTM)) . (|T97l) and (|T98| . given the relations ([210]) and (|2TT|) . respectively are 



We have set for simplicity /'(0) 



= 1, otherwise we have to renormalize the coupling constant X in the action (1 1 8 D - 
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.9 t ( t 2) (i,x) 



3 t 1 |x| 2 e- m «(l +mg) sinhm|x| 
2£ + m 2 ^ ~ 2£3 



1 , JXQe- 
W 3 

,( 3 ) 



i|x| 



TO 2 £ 3 

e(|x|-o 



=0 



(214) 
(215) 



Sift*)*) = - r 9 



(1 + m£)e 



— m£ 



2£ 3m 2 £ 3 2£ 3 



3m 2 £ 3 



2F(x) 



sinh ? 



F(0( 



1 



c| 3m 2 |x| 3 3 \|x| m|x| 2 m 2 |x| 3 
2(1 + m£)e~ m £ ( sinhmjx 



-m|x| 



e«-|x|) 



e(|x|-e)U 



to 2 £ 3 \ m|x| 

_2 2F(Q / 1 _J3 

2 |x| 3 3 V|x| m|> 



F(x; 



©(e-N) 

3 



— m|x| 



©(|x|-0 



(216) 



We note that the corrections to GR behavior are ruled by 5(x, x'). If we perform a Taylor expansion for m|x| <C 1, 
we have 



sinhi 



m x 



+ 



(217) 



(2) (2) 

For fixed values of the distance |x|, the solutions g\ t and g\j depend on the value of the radius £, then the Gauss 
theorem does not hold also if the Bianchi identities hold, as already said above [12j |. In other words, since the Green 
function does not scale as the inverse of distance but has also an exponential behavior, the Gauss theorem is not 
verified. We can affirm: the potential does not depend only on the total mass but also on the mass-distribution in the 
space. We can write 



lim F(mt;) = oo 



(218) 



obviously the limit of £ has to be interpreted up to the maximal value where the generic position |x| in the space is 
fixed. If we consider the limit £ — > (the point-like source limit), we obtain 



lim F(m£) = 1 



(219) 



By introducing three metric potentials ( I > (x), ^(x) and A(x) (the dimension is the inverse of length) we can rewrite 
(gSll and flgHj} as follows 



St ( t 2) (*,x) = r 3 $(x) 



4 2) (*,x) = r 9 *(x)% + r g A(x)^ 



(220) 



and with a fourth function, S(x), (the dimension is the cubic inverse of length) the Ricci scalar (|212j) is 



RW(t,x) = r 9 S(x) 



(221) 
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The spatial behavior (I212[) is shown in FIG.[TJ The metric potentials are shown in FIGs. [2] [3] and |U It is interesting 
to note as the function $ assumes smaller value of its equivalent in GR, then in terms of gravitational attraction we 
have a potential well more deep. A such scheme can be interpretable or assuming a variation of the gravitational 
constant G or requiring that there is a central greater mass. These two affirmations are compatible on the one hand 
with the tensor-scalar theories (in the which we have a scaling of gravitational constant) and on the other hand with 
the theory of GR plus the hypothesis of the existence of the dark matter. In particular, if the mass distribution takes 
a bigger volume, the potential increases and vice versa. 

o.oo 

-0.02 
-0.04 
£ ] -0.06 
-0.08 
-0.10 
-0.12 



2 4 6 8 10 

FIG. 1: Plot of dimensionless function £ 4 m _3 E for £ = m£ = 0.5 representing the spatial behavior of Ricci scalar at second 
order. In GR we would have S(x) = p-6>(£ — |x|). 
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2 4 6 8 10 



W/f 

FIG. 3: Plot of metric potential £m _ll 3/ vs distance from central mass with £ = m£ = 0.5. The dashed line is the GR behavior 
(similar to metric potential 




2 4 6 8 10 



M/f 

FIG. 4: Plot of metric potential £m -1 A vs distance from central mass with £ = m£ = 0.5. In GR such a behavior is missing. 

An important remark has to be done at this point. Now we can check the compatibility of /(i?)-gravity with respect 
to GR in the limit (f(R) —> R). Since the Gauss theorem is not verified for /(i?)-gravity, while the relations 
satisfy it, we have to consider the relations (I222[) and not (|214l) . (|216[) . After making the limit f(R) — > R we have 

' i?( 2 >(t,x) = 

< g\f{t^) = -ft (223) 
gg\t,x) = 

which suggest that the /(i?)-gravity is compatible with respect to GR. It is interesting to note that also in the case 
of extended spherically symmetric distribution of matter, when we perform the limit f(R) —> R, the solutions (|214|) 
and (|216p directly converge (in the vacuum) to solutions (|223[1 . showing the validity of Gauss theorem in GR. 

Another important consideration is about the asymptotic behavior of /(i?)-gravity with respect to GR. In fact, 
increasing the distance from the central mass, the gravitational field should converge to that of GR. Such a convergence 
is the standard consequence of the spherically symmetry of the source with asymptotically flat boundary conditions. 
In FIG.Q] we report the spatial behavior of Ricci scalar (|212l) approximating asymptotically the given value in GR. In 
fact supposing f(R)- gravity, the Ricci scalar acquires dynamics, and in the Newtonian limit, we find a characteristic 
scale length (m _1 ) related to the scalar massive mode. Only for distances larger than m -1 , we recover the outcome 
of GR, that is R = 0. The metric potentials are shown in FIGs. 2, 3 and 4. 

To conclude this section, we show in FIG. the comparison between gravitational forces induced in GR and in 
/(i?)-gravity considered in the Newtonian limit. Obviously also considering forces, we could obtained an intensity 
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FIG. 5: Comparison between gravitational forces induced by GR and /(i?)-gravity with £ = m£ = 0.5. The dashed line is 
the GR behavior. 



2. The Newtonian Limit of f(R) -gravity in oscillating regime 
If we consider m 2 < (i.e. sign[/'(0)] = sign[/"(0)]) from (|204j) we can choose the "oscillating" Green function 



e(x )X ') 



1 costo|x — x'| + sinm|x — x'| 

47T |x - X'l 



(224) 



The Ricci scalar (|194j) and the it-component of at 0(2) order (|196j) become 



smm x 



0« - |x|) - 2 r g ^G(0 C ° Sm|X + , Smm|X| 6(|x| - o 



(225) 
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2i?(£) sinm|x| 


m 2 £ 3 
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m 2£3 TO 


x| 




cos m|x 


+ sinm|x| 


9(| 


3 




|x| 



6(C-|x|) 



(226) 



where we introduced two new shape functions 



G(x) = 3 



mx cos mx — sm mx 



H(x) = (1 — mx) cosmx + (1 + mx) sin ma; 



(227) 



with the properties lim^o = —1 and limj _>o = !• Since we have an oscillating Green function which 

is not asymptotically zero, the "gravitational potentials" (|226|) at infinity are zero a part a possible constant value 
(lim a _>oo 2r g mG(£)(sin ma — cos mo)). 

The spatial behavior of Ricci scalar (|225[) and metric component (|226p are shown in FIGs. [5] and [7] The previous 
considerations hold also for the solutions (I225[) - (|226[) . The only difference is that now we have oscillating behaviors 
instead of exponential behaviors. 

The correction term to the Newtonian potential in the external solution can be interpreted as the Fourier transform 
of the matter density p(x). In fact, we have 



d 3 x' 

(27T)* 



p(x')e- 



-ik-x 



MG(\k\Q 
(2tt) 3 



(228) 



45 



0.4 
0.3 

M 

> 0.2 
0.1 
0.0 



10 20 30 40 

M/f 

FIG. 6: Plot of dimensionless function £ 4 m _3 E with £ = m£ = .5 representing the spatial behavior of Ricci scalar at second 
order in the oscillating case. 
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FIG. 7: Plot of metric potential ^m _1 $ vs distance from central mass with the choice £ = m£ — 0.5 in the oscillating case. 
The dashed line is the GR behavior. 

and in the limit of point-like source 

/W 3 y' , M 

-^p(x')e- lkx = (229) 

Also in this case we show in FIG.[8]the comparison between gravitational forces induced in GR and in /(i?)-gravity 
in the Newtonian limit. Obviously also in this last case we obtained a different force with respect to GR. 

3. Solutions at O (4)-order 

The metric potentials and the function S(x), respectively defined in (|220p and (I22ip . satisfy the following properties 
with respect to derivative of coordinate Z-th 9 in the matter 





We remember that |x| ; = |x| xi- 
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FIG. 8: Comparison between gravitational forces induced by GR and /(ii)-gravity with £ = m£ = 0.5 in the oscillating case. 
The dashed line is the GR behavior. 



S ; (x) 



^+(i|pl) e -™«F(x) 



Z; = $ (x)x; 



1 (l+m£) — m£ (m 2 |x| 2 +6) sinh m|x|+m|x| (m 2 |x| 2 — 6) cosh m|x| 
__ ... _ g - - - 



m 5 |x| 5 



xi = il! (-x)xi 



(230) 



A, fx) = 2 ( 1 + m ?) c -m£ (4m 2 |x| 2 +9)sinhm|x|-m|x|(m 2 |x| 2 +9)cosh m |x| = . / 
, t V, / f3 m 5 |x| 5 ' ^ 



S, In (x) = Bo(x)g, w + Ml+pil ( m 2 |x| 2 + 3) S in h: n|x|- 3m |x|cosh m |x| ^ = 3^)^+3^)^ 



_ $i„(x) = $ (x)(5 to 
and in the vacuum 



£g&e-"* ("''Ixl'+^sinhmlg-amlxlcoshmlxl ^ = $ (x)^ + $1 ( x )x,X n 



A,;(x) 



|x| 5'^ F(Qe- m W Xl = S (x)x ; 

_J_ m|x| + lF(Q e — = $ , x 
| x |3 T | x |3 3 x t — Vol^Xj X; 



n 2 |x| 2 -2 _ m J |x| J -m 2 |x| 2 -6m|x|-6 F(£) e - m|x| 
m 2 |x| 5 m 2 |x| 5 3 



m 2 |x| 5 



m J |x| J +4m 2 |x| 2 +9m|x|+9 2F(g)e~ 
m 2 |x| 5 3 



Xi = * (x)x; 

.T; = A (x)xi 



(231) 



S, In (x) - Bo(x)g, w - r^S^^^ F ^) e -m\ x \ XlXn = So (x)<5 in + S 1 (x)x i x n 



$i„(x) = $ (x)(5 /n 



3 , m 2 |x| 2 +3m|x|+3 F(g)e~ 



TxP 



X;X„ = $o(x)<5i„ + $i(x)x;X„ 



Obviously when we consider the physics in the matter or in the vacuum we have to choose the "right" quantities 
Ho(x), Hi(x), $ (x), $i(x), * (x), Ao(x). 
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The expression of Ricci scalar at fourth order (|199[) is 



2m x 



dlx'llx'M e" m||x| ^ |x '" - e - m (l x l+l x 'D 



,'|2 



3A(x0 + W^O±A« |xT ' 



So(x') + m 2 *(x')S(x') + S 1 (x')A(x')|x'| 2 (232) 



where we note two contributions. The first one still depends on the quadratic term (oc R 2 ) in the action (|18[) . while 
the second one is related to the cubic term (oc R 3 ). By introducing the two functions Sr(x) and S/j(x), Eq. (|232[) 
can be rewritten as follows 



i—i / S i—i / \ 



(233) 



An analogous situation is found for the tt-component of metric tensor at fourth order. In fact Eq. (|201[) becomes 
r„X /* 



3t ( t 4) (t,x)= d|x'||x'||||x|-|x'||-|x|-|x'||p(x')$(x') 
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2ro 2 $(x') + *(x') -3$ (x')-*i|x 
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M |X| Jo 



and by introducing other new functions $/(x), <I>jj(x), finally we have 



3m 4 



A x ,5 7J (x') 
3m 2 



(234) 



9t ( t %x)=r g 2 



*iW + -r*//(x) 



(235) 



It is useful to note that we have generally four contributions to g^r in (I234[) . The first one is induced by the 
non-linearity of the metric tensor even in static spherically symmetric case. The product p(x)<I>(x) is non-zero only 
in the matter but contributes to determination of tt-component in any point of the space. The second one takes into 
account the contribution induced by the solution of previous order for the determination of the ii-component of the 
Ricci tensor at fourth order. These first two terms are present also in GR. While the second two terms are derived 
from the modification of the theory. In fact the third contribution depends on the addition of the quadratic term 
(oc R 2 ) in the action and finally the fourth one from the addition of the cubic term (oc R 3 ). 

The choice of free parameter /z, which is linked to the third derivative of f(R), is a crucial point in both the 
expressions (|233|) and (|235[) to obtain the right behavior. From the mathematical interpretation of Newtonian limit 
one has |/"'(0)| < |/"(0)| and if p 4 > (i.e. sign[/'(0)] = -sign[/"'(0)], otherwise [i 4 is not a length) we have 
m 4 //i 4 — |/"'(0)|/3/"(0) 2 , so we find the constraint < m 4 /fi 4 < 1. In FIG. we report the spatial behavior 
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of (|233l) in the matter and in the vacuum, (0.3 ^ m 4 /fi 4 ^ 0.9), showing that far from the source we obtain a 
spacetime with a vanishing scalar curvature. At Newtonian level, the Ricci scalar (|212l) is negative defined while 
at post-Newtonian limit it is positive defined. 

In FIG. (TTU1) . we report the time-time component of metric tensor g[f* on the same interval of values of m 4 //i 4 , 
although the behavior is quite insensitive to changes induced by the contributions of the cubic term in the Lagrangian. 
Besides we can observe an important analogy with respect the results of GR. In both cases, we have a potential barrier, 
but for /(i?)-gravity it is higher (as in the Newtonian limit we found a deeper potential well). 




2 4 6 8 10 

W/f 



FIG. 9: Plot of dimensionless function £ 4 m 4 (Hj + z£ii), representing the Ricci scalar at fourth order, where z — m 4 /^ 4 
and £ = m£ = 0.5. The spatial behavior is shown for 0.3 ^ z ^ 0.9 (solid lines) while the dotted line corresponds to 
R - tt^R 2 -theory. 




FIG. 10: Plot of dimensionless function <^ 2 m _2 ($/ + z&u) (solid lines) where z = m 4 /fi 4 and of function l/2|x| 2 (dashed 
line). For z = (dotted line) we have the behavior of R — g^-ii^-theory. The solid lines are obtained for 0.3 ^ z ^ 0.9 and 
for £ = m£ = 0.5. 



4- Solutions from isotropic to standard coordinates 

The metric solutions that we have found are expressed in isotropic coordinates and often, for spherically sym- 
metric problems, they are conveniently rewritten in standard coordinates (the standard form in which we write the 
Schwarzschild solution). Here the relativistic invariant of metric (|189[) is 



ds 2 



l + r 9 $(x) 



$/(x) + — $jj(x) 



dt 2 - 



1 - r 9 *(x) 



|dx| 2 + r 9 A(x) 



(x • dx) 



(236) 
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From spherically symmetric form of (|236l) . it is convenient to replace it with Eq. (|42|) by transformation (|34|) . The 
proper time interval (|236[) then becomes 



ds 2 



1 + r g $(r) + r 2 $ z (r) + -^-$/r(r) 



dt- 



l-rj *(r) + A(r) 



dr 2 



1 - r g *(r) 



r 2 dn (237) 



To get the metric in the standard form, we need to impose a radial coordinate transformation 



1 - r g ^(r) 



2 ~2 

r = r 



(238) 



and we have a new set of coordinates {r, 9, 0}. The metric (|237[) becomes 



ds 2 = 



1 + r g $(f) + r g z ( $/(f) + — $//(f) 

A* 



dt 2 - 



1-rJ *(f) + A(f) 



^ ) df 2 - f 2 dfi 



(239) 



The explicit expression of (|239[) is not displayed because the equation (|238l) cannot be solved algebraically. However, 
this technical problem is overcome with the help of numerical methods when we are interesting to test experimentally 
the theory. 



VIII. THE NEWTONIAN LIMIT OF QUADRATIC GRAVITY: f(X,Y) = aiR + a 2 R 2 + a s R aP R al3 

Since terms resulting from R n with n > 3 do not contribute in the Newtonian limit, as we have seen previously, 
we provide explicit solutions for different types of Lagrangians generated by a primitive Lagrangian, the so-called 
Quadratic Lagrangian of the form [58} 



f(X, Y) = a x R + a 2 R 2 + a 3 R a[j R a ^ (240) 

where oi, a 2 and 03 are arbitrary constants 10 . Such a Lagrangian belongs to the general class of FOG (fTE)) . The field 
equations coming from (|240p arc 

H^v = (ai + 2a 2 R)R llv — aiR+a2R + a 3R a pR — _ 2a 2 R tfll , + 2a 2 g^DR + 2a 3 R^ a R au — 2a 3 i? Q ( /i;! ,) Q 

+a 3 DR fiU + a 3 R a ^g„ v = X (241) 

, H = -axR + 2(3a 2 + a 3 )DR = XT 

Till now the solutions of field equations (in general Eqs. (|2"I]>) are found by considering the trace Eq. (j2"2")l as the 
dynamical equation for the Ricci scalar. This approach allows to study a set of second order differential equations 
while the FOG equations are intrinsically fourth order differential equations in metric formalism. By using the trace 
equation (second line of (|241jl ) or the definition of the Ricci scalar ([5]) and substituting them into the field equations 
(first line of (|241jl ). we have a set of fourth order partial differential equations. In this section, we will study FOG 
field equations. 

If we introduce the gravitational potentials in the isotropic metric (14"TT) by the quantities $ and $ linked to g^' 
and gry 



ds 2 = 



1 + 2$ 



dt 2 - 



1 - 2* 



5ijdx l dx 3 



(242) 



Note that the physical dimensions of the constants are [02] = [123] = length 2 and [fti] = length . 
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and by using the paradigm of Newtonian and post-Newtonian developments, we can investigate, without assuming 
the harmonic gauge condition, field Eqs. (|24ip for the Quadratic Lagrangian (|240[) at Newtonian order, that is 



2aiA* - 2(4a 2 + a 3 )A 2 * + 2(2a 2 + a 3 )A 2 $ = Xp 



A 



ai(* - $) + (4a 2 + a 3 )A$ - (8a 2 + 3a 3 )A* 



ai(* - $) + (4a 2 + a 3 )A$ - (8a 2 + 3a 3 )A* 



(243) 




By introducing two new auxiliary functions ($ and Eqs. (|243|) become 



2,1 2 



2a 2 +a 3 



A 2§ _ 2 03 (3a 3 +a 3 ) A 2^ _ 4a2±a3 A< | _ «^ A ^ = ^ p 
ai(2a 2 +a 3 ) 2a 2 +a 3 2a 2 +a 3 r 



A 



$ + A* 



$ + A* 



(244) 



where $ and $ are linked to $ and $ via the transformations 



3a 2 +3a 3 )$+Qi1' 
2ai (2a 2 +a 3 ) 



2ai(2a 2 +a 3 ) 



(245) 



Obviously we must require a!(2a 2 + a 3 ) 7^ 0, which is the determinant of the transformations (I245p . Let us introduce 
the new function 3 defined as follows 



n=$+ a* 



At this point, we can use the new function S to uncouple the system (I243p . With the choice $ 
possible to rewrite Eqs. (|243[) as follows 



(246) 
3 - A*, it is 



2a 3 (3a 2 +a 3 ) A3 ^ 602+03 A 2§ _ aj_ A y = X p + T 



AJtfy-^y = 



(247) 



where r = 4a 2 +a 3 A j + 2a 3 (3a 2 +a 3 ) Wg are interested in the solution of (EH in terms of the Green function 

2a 2 +a 3 ai(2a 2 +a 3 ) Y ' 

~l(x,x') of field operator ^^p|^A 3 + ff|±f|A 2 - ^fe A - Thcn Eqs.flMS]) are equivalent to 



( 20^302+03) A 3-, ( m + 602^0, A 2-,/ A _ oj_ A -[( = <J( X _ X ') 

ai(2a 2 +a 3 ) V / 2a 2 +a 3 'V ' / 2a 2 +a 3 \ ' / V J 



(248) 



[ A3(x)^-3(x)^ = 



The general solutions of Eqs. (|244l) for $(x) and ^(x), in terms of the Green function ~I(x, x') and the function 
J(x), are 



$(x) 



So 2 +3a 3 )A x — ai r j3„/ 



2ai(2a 2 +a 3 



Jd 3 x'l(x,x') 



A>(x') + |^±^A X ^(X') + 2a 3(3a 2 +a 3 ) A 2 ^ 
r\ I 2a 2 +a 3 x \ > ai(2a 2 +a 3 ) x""V I 



8a 2 +3a 3 
2ai(2a 2 +a 3 ) 



J(x) 



(249) 



402 +03 



2ai(2a 2 +a 3 ) 



Eqs. (|244[) represent a coupled set of fourth order differential equations. The total number of integration constants 
is eight. With the substitution (|246l) . it has been possible to decouple the set of equations, but now the differential 
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order of the single equations is changed. The total differential order is the same, indeed we have an equation of sixth 
order and another equation of second order, while, previously, we had two equations of fourth order. The number 
of integration constants is preserved. We can conclude that, with our approach, also the introduction of the new 
quantities $, "J does not contradict the paradigm of metric theories of FOG. The price is that now the r.h.s. of 
tt- component of the field equation has been modified: there is an additional matter term r coming from the ij- 
componcnt. In Table HTT| we show particular cases of (|244|) for different choices of coupling constants of the theory 
with the vanishing determinant of transformations (|245[) . 



Case 


Choices of constants 


Corresponding field equations 


A 


a 2 = 
a :i = 


A* = ^p 

a[$(x)+ g frf 3 xV (x '!, 

V / 1 a\ J |x — x' 


A JO - — 


$(x)+ G fd 3 x', p{x ']] -0 

J ,i.7 


B 


ai = 
a 3 = 


A 5 
A 


(2* - $) = 
A(2* - $) 




A(2* - $) 


= 


C 


ai = 
a 2 = 


A 5 
A 


(*-*) = 2 
A($ - 3*) 


(Jo T — 


A($ - 3*) 


= 


D 


a3 = — 2a 2 


2a 2 A 2 * - ajA* = -fp 

V 2 ai$(x) 2a 2 A$(x) + G J d 3 x' 




ai$(x) 2a 2 A$(x) + G/d 3 x , | p(x '! | -0 

J ,17 


E 


oi = 
a 3 = -4a 2 


A : 
A 


* = 
A* 


A 1 „ 

- 


A* 


= 


F 


ai = 
03 = — 2a 2 


A A(* - $) 




A(* - $) 


= 


G 


oi = 


A : 
A 


(2* 
A$ 


+ $) = 

5y - 


3. 
4c 

A$ 


= 

Ai 



TABLE III: Explicit form of the field equations for different choices of the coupling constants for which the determinant of 
transformations (|245l) vanishes. 



A. Green functions for systems with spherical symmetry 

As above, we are interested in the solutions of (|244[) at Newtonian order by using the method of Green functions 
with spherical symmetry. Let us introduce the radial coordinate r == |x — x'|; with this choice, the first equation of 
(|248|) for r ^ becomes 



2a 3 (3a 2 + a 3 )A 3 r 1(r) + (6a 2 + a 3 )A 2 r l{r) - a\ A r 1(r) 







(250) 



where A r . 



2 d r (r 2 d r ) is the radial component of the Laplacian in polar coordinates. The solution of (|250l) is 



l{r) = Kx - - 

r 



K 2 • —I 

ai 



2(3a 2 + a 3 ) 



oi 



V 



2(3a 2 + a 3 ) ' 



(251) 
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where Ki are constants. We note that, if a 2 = 03 = 0, the Green function of the Newtonian Mechanics is recovered. 
It is the same situation of the Electromagnetism. The integration constants Ki have to be fixed by imposing the 
boundary conditions at infinity and at the origin. In fact Eq. (|251[) is a solution of (|250|) and not of the first equation 
in (|248[) . A physically consistent solution has to satisfy the condition ~I(x, x') — > 0, if |x — x'| — > 00, then the 
constants K\, K4, Kq in Eq. (125ip have to vanish. To obtain the conditions on the constants K2, K3, K$, we consider 
the Fourier transformation of ~l(x,x'), that is 



"x.x'j = /^1(k)e*("'' (252) 



By putting (|252[) in the first equation of (|248[) . we obtain 



j3u ik(x-x') 

~l(x x ) = - / - (253) 

/ (27r) 3 2a 3(3a 2 +a3h t 6 _ 602+03 ^4 _ 01 J £ 2 
01(202+03) 2a 2 +a 3 2a 2 +a 3 



which, in the case of spherical symmetry, becomes 1 



1(x,x') 



1 ai(2a 2 + a 3 ) 1 
An 2 a 3 (3a 2 + a 3 ) |x — x' 



d|k|sin|k||x- 



k 2 _ 21 



k 2 



2(3a 2 +a 3 ) 



(254) 



The analytic expression of ~I(x, x') depends on the nature of the poles of |k| and on the values of the arbitrary 
constants ax, <x 2 , a 3 . If we define two "masses" m\ and m 2 



ax 



2(3a 2 + a 3 ) ' 



ax 

«3 



(255) 



we obtain, in Table IIV[ the three particular expressions of (|254|) . We note that the definition of mass mi is a 
generalization of (|195[) (and obviously of (|154[) and (I185[) ^) . 



Case Choices of constants 



A 



£13 > 

3a2 + a3 < 



a 3 < 

3a2 + £13 > 



£13 < 



3(J2 + £13 < 



Green function 



n A (x,x') 



-|B( X)X ') = -rJ-r-^TT 

' / YZtt Ix — x 



cos(mi |x — x'|) 



s(m 2 \x— x^l) 



n c (x,x') 



I cos(m2 Ix— x' |) 

H 



TABLE IV: The complete set of Green functions for (|254[1 . It is possible to have a further choice for the scale lengths depending 
on the other two length scales. If we perform the substitution m\ ^ — mjj, we obtain a fourth choice. In addition, for a correct 
Newtonian component, we assumed a\ > 0. In fact when 02 = £13 = the field Eqs. (|243p give the Newtonian limit for 
ai = 1. 



Besides, we have obtained another Yukawa-like correction to the Newtonian potential related to the squared Ricci 
tensor correction in the Lagrangian (I240|) . This behavior is strictly linked to the sixth order of (I248|) . which depends 



11 we introduced the polar coordinates in the k-spacc. 
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on the coupled form of the system of equations (|243l) . In fact, if we consider the Fourier transform of the potentials 
$ and * 



$(x) = 



d 3 k 

(2tt) 3 / 5 



$(k) e 4k 



*(x) = 



rf 3 k 

(2tt) 3 / 5 



*(k) e* 



(256) 



the solutions are 



$(x) 



*(x) 



2 J (2tt 



ai + (8a2+3a 3 )k 2 ]p(k)e' k - 3c 

(2tt) 3 / 2 k a (ai^a 3 k^)[a 1 +2(3Q2+a3)k^ 

* f rf 3 k [ai + (4a 2 +a3)k 2 ]p(k)e ik - 3c 
" 2 J (2tt) 3 / 2 k 2 (ai-a 3 k 2 )[ai+2(3a 2 +a 3 )k 2 ] 



(257) 



where p(k) is the Fourier transform of the matter density. We can see that the solutions have the same poles as 
Eq.flSMD. 
([2"ST|) are 



Eq. (|254[) . Finally, considering the Fourier transform of the point-like source (fSTj) , that is p(k) = p^Z™ , the solutions 



$( X ) = — GM 1 1 eZ^W _ 4e^ 
V / ai I |x| 3 |x| 3 |x| 



GM I 1 1 e ~ mi|3c| _ 2 e"" 12 
3 |x| 3 — [xf 



(258) 



where 03 ^ and 3a2 + 03 7^ and mf , m| > 0. The solution for $ generalizes the third line of (|222l) in the case 
of point-like source. In fact, it is interesting to note that, if 03 = (m2 — > 00), we have the missing of a scale length 
(a pole is missed) with only a Yukawa-like term analogously to Electrodynamics. The Green function, in this case, is 



1 < k »— " wCT (25<J) 

and Lagrangian (|240l) becomes / = a\R + a 2 R 2 . Finally the presence of the pole is achieved considering a particular 
choice of the constants in the theory: 03 = — 2a2- In Table UTT1 ( Case D), we provide the field equations for this choice 
and the relative Green function is 



W^Wkjoc ^^^ (260) 

The spatial behavior of (|259[) - (|260[) is the same but the coefficients are different since the theories are different. In 
conclusion, we need the Green function for the differential operator A 2 . The only possible physical choice for the 
squared Laplacian is 



V)(x.- X ')cx^— ^ (261) 

since the other choice is proportional to |x — x'| and cannot to be accepted. Considering the last possibility, we will 
end up with a force law increasing with distance [17] . In summary, we have shown the general approach to find out 
solutions of the field equations by using the Green functions. In particular, the vacuum solutions with point-like 
source have been used to find out directly the potentials, however it remains the important issue to find out solutions 
for systems with extended matter distributions. 



B. Solutions using the Green functions approach 



Before investigating the general solution of (|243|) . we want discuss all cases shown in the Table UTTl Later, we will 
derive solutions in presence of matter using the Green functions in Table IIVI 
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Specifically, in Table [Vj we provide solutions in terms of Green functions of the corresponding differential operators 
coming from the field equations shown in Table IIIIl Case A corresponds to the Newtonian theory and the arbitrary 
constant a\ can be absorbed in the definition of matter Lagrangian. The solutions are 



i$(x) = A *(x) 



G / <f x 



3,./ />( x ') 



(262) 



For Case D, instead, we have the field equations of a sort of modified clcctrodynamic-like representation. The solution 
can be expressed as follows 



D*(x) = D *(X) = -G / dV 



rl*-*'h 



(263) 



The solutions make sense only if ai/ag > 0, then we can introduce a new scale-length. A particular expression of 
(|263|) . for a fixed matter density p(x), will be found in a more general context in the next section. Nevertheless these 
two cases are the only ones which exhibit the Newtonian limit (obviously the first one!), while for the remaining cases 
there are serious problems with divergences and incompatibilities. In fact, Case B presents an incompatibility between 
the solution obtained from the it-component and the one from the ^/'-component. The incompatibility can be removed 
if we consider, as the Green function for the differential operator V 4 , the trivial solution: G(a 2 )\b — const. With this 
choice, the arbitrary integration constant <J>o can be interpreted as —GM. However another problem remains: namely 
the divergence at the origin. The interpretation of the constant $o as a total mass and not as a generic integral 
J cZ 3 x'p(x') does not avoid the singularity. Finally, the solution 



2 B *(x) - B $(x) 



GM 



(264) 



holds only in vacuum. Before continuing the analysis of the various cases, the term J d 3 x'~I(^2)(x, x')/o(x') has to be 
discussed for the choice (|261l) . A generic field equation with A 2 (from Table UTi]) is 



A 2 $(x) ex A 2 J d 3 ^'-0^- « A xP (x) + 47rp(x) 



(265) 



from which the only consistent possibility is to set p(x) = 0. In the remaining cases, we can only consider vacuum 
solutions. 

By considering solutions (|249j) with the Green function "I (x, x') from Table |IV] and by assuming 3(x) = 12 , we 
have 



(266) 



After some algebraic calculations we get, for a spherically symmetric matter source (I203p . the internal and external 
solutions for $ 



t$m(x) 



l$o«t(x) 



GM 
ai 



GM 



8mf +m|(2+3mfg 2 ) 



J_ F(miQ e 
|x| + 3 



J5l! 
2?3 



(l+mi^) sinhmi|x| 
mi|x| 



+ 4- 



|x| 



4F(m 2 C) e - 



|x| 



' 12 ^tl+ m 2C) sinhm2|x| 
m|£ 3 m 2 |x| 



(267) 



12 We have to note that our working hypothesis, 3(x) = 0, is not particular, since when we considered the Hilbert-Einstein Lagrangian to 
give the Newtonian solution, we imposed an analogous condition. 
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Case 


Solutions 


Newtonian behavior 


A 


$(x) - *(x) - ° f d 3 x', p(x '' 

\ / V ' a\ J x — x'| 


yes 


B 


2*(x) - *(x) = 

2*(x) - $(x) = Jd 3 xH (V 4 ) (x,x')p(x') 


no 


C 


*(x)=^ + ^/d 3 x'-I (v4) (x,x') /9 (x') 


no 


D 


$(x) = -4 7 rG/d 3 xH (2(l2V 4_ aiV 2 ) (x,x')p(x') 
*(x) = -4ttG/ d 3 x / 1 f20!lV 4_ 0lV2 ,(x,x')p(x / ) 


yes 


E 


$(x) = -^/d 3 x'l {V 4 ) (x,x')p(x') 
*(x) = ft 


no 


F 


$(x) = ft - f JrfVl (v4) (x,x>(x') 
>t(x)^f /d 3 x'l (v4) (x,x')p(x') 


no 


G 


$(x) = ft 

*W = "lit -f /d 3 x'V ) (x,x')p(x') 


no 



TABLE V: Solutions of the field equations in Table Hill Solutions are found by setting 3 — in the ij -component of the field 
Eqs. (|248p . Solutions are displayed in terms of the Green functions. $o is a generic integration constant. 



representing a generalization of metric potential (j2 14|) . It is easy to check that when a^, — > (i.e. mi — > oo), we get 
exactly Eq. f)2 14|) . If we consider the limit of point-like source, we get the first line of (|258j) . 

For the sake of completeness, let us derive solutions for the other Green functions. Starting from Case B in Table 
IIV1 yet we have the internal and external potentials 



**<»(*) = -SM. 
B <I> out (x) = 



x r> (2-\-3m 1 $ > )—8m 1 |x| cos mi^+mig sin sin mi |x| , a cos m 2 f;+m2£ sin m^j sin m 2 |x| 

2mfm|£ 3 2^ 3 mi|x| 4 m 2 |x| 



1 . G(mi^) cos mijxj 4G(m.2g) cosm.2|x| 
R+ 3 |x| 3 |x| 



(268) 



and its limit of point-like source is 



lim s$ „t(|x|) 



GM 



a i 



1 lcosmi|x| 2cos?7i2|x| 
1x7 ~ 6 bd + 3 ixj 



(269) 



Finally for Case C in Table ITVl we have 
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c<Mx) = 



t 2 (2+3m^ 2 )+8mj _ |x| 2 _ e "^(l+m^) sinhmi|x| _ / cos rre 2 £+TO 2 g sin m 2 C ginmgjxj 
2m^m^ 3 2£3 ^2|3 ^TJ| 1 ^Z|3 m 2 |x| 



1 

|x| 



c$ out (x) = 
which in the point-like limit becomes 



^("HC) e~'" l|x| , 4G(to 2 Q cosm 2 |x| 
3 |x| 3 |x| 



(270) 



lim c $ otlt (x) = 



GM 



1 



le -mi|x| 4costo 2 |x| 

3^d 3 bd 



(271) 



Similar relations are found for ^ . 

These results mean that, for suitable distance scales, the validity of Gauss theorem is restored and the theory agrees 
with the standard Newtonian limit of GR. 



IX. THE NEWTONIAN LIMIT OF f(X, Y, Z)-GRAVITY 



We conclude this review providing the Newtonian limit of a generic f(X, Y, X)-gravity [5£ 
paradigm of Newtonian limit (1981) . the curvature invariants X, Y , Z become 



By considering the 



r x ~ x^ + xw 



(272) 



and the function / can be developed as 

f(X,Y,Z) ~ /(0) + /x(0)X( 2 ) + i/ XX (0)X( 2 ) 2 + / x (0)x( 4 )+/ y (0)r( 4 )+/z(0)Z( 4 ) + ... (273) 
and analogous relations for partial derivatives of / are obtained. From the lowest order of field Eqs. ([2~Tj) . we have 



/(0) = 



(274) 



this means that not only in /(i?)-gravity [36|, [45] but also in f(X, Y, Z)-gravity, a missing cosmological component in 
the action (fl~8f implies that the space-time is asymptotically Minkowskian. Eqs. (f2~Tj) and ([22]) at C(2)-order become 13 



= /x(0)4 2) - [/y(0) + 4/z(0)]A4 2) - M^X( 2 ) - [/ XJC (0) + ^]AI< 2 ) = #T t ( t 0) 

= f x (0)R^ - [/y(0) + 4/ z (0)]Ai? 4 (2) + MX( 2 )% + [/ XJf (0) + IzM]AXW5 l3 - / XX (0)X( 2 ),^ + 



+[/r(0) + Af z (0)]R% m + f Y (0)R%{ im = 
# (2) = -fx(0)XW - [3f xx (Q) + 2jy(0) + 2/ z (0)]AX( 2 ) = XT^ 



By introducing the quantities 



(275) 



13 We used the properties R a p> ari = §QR and R^ v . a n = R^ a . ua ~ OR^ 
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m 2 j_ fxM 

" ll ~ 3/xx(0)+2/y(0)+2/ z (0) 



m 2 2 = 



/x(0) 
/y(0)+4/ z (0)] 



we get three differential equations for curvature invariant X^ 2 \ tt- and ij-component of Ricci tensor R 



(2) 



(276) 



(A - m 2 2 )R$ 



\ (A-m 2 2 )R < ? ) + 



r»2 _ mi +2m 2 A 
2 



m 2 2 -y t(0) 

"/x(0) J « 



mi — rn 2 f)2 



12_ _ mi +2m 2 A W 
2 6ml- 1 J L J 



xm = o 



(A-mi 2 )X (2 » = f^yTC°) 



(277) 



We note that the definitions (|276[) are a generalization of (|255p . While the interpretation of rci\ is clear from the trace 
field equation of /(-R)-gravity, now also m 2 is clear. In fact, from the two first lines of (|277j) the Ricci tensor has a 
characteristic length (?ti 2 _1 ). 

The solution for curvature invariant X^ = R^ at third line of (|277[) is 



*( 2 >(t,x) = ™/ d 3 x^ 1 (x,x')T(°)(t,x') 



(278) 



/x(0) 

where 0i(x, x') is the Green function of the field operator A — mi 2 (this solution is formally equal to (|194[1 ). The 



solution for g^' , by (|103[) . is the following 



m?X (Q) „ (m 1 2 + 2m 2 2 )X (0) m 2 2 - m x 2 m 

x(0) " 1 ' j 6/ x (0) ' )+ 6 A j 

(2) 



279) 



where C/ 2 (x, x') is the Green function of field operator A — m 2 2 . The general solution for from Eqs. (|277[) . if we 
consider the expression (| 1 16|) (gauge harmonic condition), is 



mi 2 - to 2 2 2 / to 2 2 mi 2 + 2m 2 . . . 

O.v/i// - — : s A x // d. 



3mi 2 V 2 6mi 2 ^ x ' '' 

?(2) 



X^(x.") (280) 



If we choose the metric (|242p . from (I103[) we have R^ = A^d^j + — 3?),-y and the second field equation of 
(12TT)) becomes 



A* = /We 2 (x,x') ^^B^A X , I< 2 )(x') 



(281) 



(*-*),« = ^ r %^/d 3 x'&(x,x')X^,(x') 

(2) ______ _____ 

Then the general solution for ^ from (|277p . without gauge conditions and by using the first line of (|281D . is 



si? =2^Sii 



2ir 



7 |x-x'| V 2 6mi 2 y v ; 



2 6m 

and the second line of (|281l) is only a constraint condition for metric potentials. In fact from its trace, we have 



(282) 



A(*-*) = I '/ :i x'c7.,(x.x'l^ x .A^'(x'l 



(283) 
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and we can affirm that only in GR the metric potentials $ and ^ are equals. 

Let us consider the point-like source (jSTj) . If we choose mi 2 > and m 2 2 > 0, the curvature invariant X^ (|278[) 
and the metric potentials $ (|279|) and * (|282|) are 



$ = 



" fx(0) |x| 



GM 
7x(0) 



* = - 



GM 

7IW) 



i 

|x| 



1 

|x| 



1£ 
3 



3 



|x| 



(284) 



1 e~ m i 
3 |x| 



2 e - m 2 

3 |x| 



The modified gravitational potential of /(i?)-gravity is further modified by the presence of functions of R a pR a ^ (as 
in the Quadratic Lagrangian-theory) and R a j3 7 sR al3jS ■ The curvature invariant X^> presents a massive propagation 
and when f(X,Y,Z) ->■ f(R) we find the mass definition m 2 = -f'(R = 0)/3/"(i? = 0) given in ([TTJSj) . In this 
case, the propagation mode with disappears. Obviously the two expressions for gravitational potentials in (|284[) 
satisfy the constraint condition (|283[) . In FIGs. [TT1and[T2lwe sketch the spatial behavior of metric potentials for some 
values of parameters m\ and m^. 




W 



FIG. 11: Plot of metric potential $ in Eqs. (|284[) . ni2 = C m i an d mi = 0.1 (dotted line), mi = (m2 and m2 = 0.1 (dashed 
line). The behavior of GR is shown by the solid line. The dimensionless quantity £ runs between -i- 10. We set fx(0) ~ 1. 




FIG. 12: Plot of metric potential ^ in Eqs. (|284[) . 7712 = C m i an< l m i = 0.1 (dotted line), mi = (m2 and mi = 0.1 (dashed 
line). The behavior of GR is shown by the solid line. The dimensionless quantity £ runs between -i- 10. We set fx(0) = 1. 

The same outcome has been found for Quadratic Lagrangian if we identify a\ — fx(0)- We can affirm, then, the 
Newtonian limit of any f{X, Y, Z)-Gravity can be reinterpreted by introducing the Quadratic Lagrangian and the 
coefficients have to satisfy the following relations 
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a 1 =f x (0), a 2 = -fxx(0)- fz(0), a 3 = fy(0) + 4f z (0) (285) 

if we want the same definition of parameters in (|255l) and (|276[) . 

A first considerations about (|285[) concerns the characteristic lengths induced by f(X,Y,Z)-theory. The second 
length 7Ti2 _1 is originated by the presence, in the Lagrangian, of squared Ricci and Riemann tensors, but also a 
theory containing only squared Ricci tensor shows the same outcome. Obviously the same is valid also with the only 
squared Riemann tensor . Such terms give rise to massive gravitational modes and then to the possibility of massive 
gravitational waves (see [6lj and references therein) . 

A second consideration concerns the Gauss - Bonnet invariant defined by the relation Gqb — X 2 — AY + Z 60] . 
In fact the induced field equations satisfy, in four dimensions, the condition 



H^ B = - + H% = (286) 

and, by substituting them at Newtonian level (Hf t ~ — AAR^) in Eqs. (|2ip. we find the field equations (ever at 
Newtonian Level) of Quadratic Lagrangian. 

A third and last consideration is about the solutions (|284[) . When we perform the limit in the origin |x| = we 
have no divergence. In fact we find 



mi — 4toq mi + 2m? 

lim $ oc — lim * oc 287 

|x|->o 3 |x|->o 3 v ; 

and only if we remove, in the action (TT5|) . the dependence on the squared Ricci or squared Riemann tensors, we get 
the divergence of GR. For a physical interpretation of solution (|284p . we must impose the condition mi — 4m,2 < 
to have a potential well with a negative minimum in |x| = and mi < m 2 to have a negative profile of potential 
(see FIG. [TTj) . Then, if we suppose fx(0) > 0, we get a constraint on the derivatives of / with respect to curvature 
invariants, that is 



fxx(0) + f Y (0)+2f z (0) < (288) 

In the case of /(i?)-gravity (/V(0) = /z(0) = 0) we reobtain the same condition between the first and second 
derivatives of f(R). 



X. CONCLUSIONS 



The weak field limit is a crucial issue that has to be addressed in any relativistic theory of gravity. It is also the 
test bed of such theories in order to compare them with the well-founded experimental results of GR, at least at Solar 
system level. 

In this review paper, we have considered the problem of weak field limit of Fourth Order Gravity, that is of 
gravitational theories where curvature invariants, a part the standard Ricci scalar, linear in the Hilbert-Einstein 
action, are taken into account. In particular, we have analyzed the Newtonian and the post-Newtonian limits of 
theories involving non-linear combinations of Ricci scalar, Ricci tensor and Riemann tensor. The calculations have 
been essentially developed in the so-called Jordan frame but we have also considered the conformal transformations 
and the possible shortcomings emerging in carrying the weak field limit results in the Einstein frame without an 
appropriate interpretation of post-Newtonian parameters. 

The general feature that emerges from the weak field limit is that corrections to the Newtonian potential naturally 
come out. These corrections are Yukawa-like terms bringing characteristic masses and lengths. Conversely, the 
standard Newtonian potential is just a feature emerging in the particular case f(R) = R. These characteristic masses 
(and lengths) come out as combinations of the parameters of the theory and fix the scales where corrections become 
relevant. 

These results open new intriguing possibilities since accurate measurements of PPN parameters could confirm or 
rule out these theories in view of the forthcoming space experiments as GAIA and GAME (see [62- 64] for a detailed 
review of the status of art). 
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On the other hand, it is well-known that the new features related to extended theories of gravity could have 
interesting applications in other fields of astrophysics as galactic dynamics (65[, large scale structure [66j and cosmology 
[H3| in order to address dark matter and dark energy issues. The fact that such "dark" structures have not been 
definitely discovered at fundamental quantum scales but operate at large astrophysical (infra-red scales) could be due 
to these corrections to the Newtonian potential which can be hardly detected at laboratory or Solar System scales. 
Finally, the presence of unavoidable light massive modes could open new opportunities also for the gravitational waves 



detection of experiments like VIRGO, LIGO and the forthcoming LISA [61 
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